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ABSTRACT 

We study the origin and properties of "extra" or "excess" central light in the surface brightness profiles of 
remnants of gas-rich mergers. By combining a large set of hydrodynamical simulations with data on observed 
mergers that span a broad range of profiles at various masses and degrees of relaxation, we show how to robustly 
separate the physically meaningful extra light - i.e. the stellar population formed in a compact central starburst 
during a gas-rich merger - from the outer profile established by violent relaxation acting on stars already present 
in the progenitor galaxies prior to the final stages of the merger. This separation is sensitive to the treatment 
of the profile, and we demonstrate that certain fitting procedures can yield physically misleading results. We 
show that our method reliably recovers the younger starburst population, and examine how the properties and 
mass of this component scale with e.g. the mass, gas content, and other aspects of the progenitors. We consider 
the time evolution of the profiles in different bands, and estimate the biases introduced by observational studies 
at different phases and wavelengths. We show that, when appropriately quantified, extra light is ubiquitous in 
both observed and simulated gas-rich merger remnants, with sufficient mass 3 - 30% of the stellar mass) 
to explain the apparent discrepancy in the maximum phase-space densities of ellipticals and their progenitor 
spirals. The nature of this central component provides a powerful new constraint on the formation histories 
of observed systems and can inform both our studies of their progenitors and our understanding of the global 
kinematics and structure of spheroids. 

Subject headings: quasars: general — galaxies: nuclei — galaxies: active — galaxies: evolution — cosmology: 
theory 



1. INTRODUCTION 

A wide range of observed phenomena support the view 
that gas-rich mergers play a central role in galaxy evo- 
lution. The most intense starbursts, ultraluminous in- 
frared galaxies (ULIRGs), are always associated with merg- 
ers (e.g. IJoseph & Wright! Il985l: iSanders & Mirabell 1199a) . 
with dense gas in their centers providing material to feed 
black hole (BH) growth and to boost the concentration and 
central phase space density of merging spiral s to match 
those of ellipticals (Hernquist et al. 1993; Ro bertson et alJ 
I2006bl) . Studies of ongoing mergers and remnants from 
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mergers of gas-ri ch systems (e.g., iL ake & Dressier] 
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Dovon et al.l ll994l; IShier & Fischer 199 
Genzel et al."2001'; Tacconi et al. 2002; iDasyra et al 
2007; Rothberg & Joseph 2004, 2006a1), as well as post- 
starburst (Eh-A/Kh-A) galaxies (Goto 2005), have shown that 
the kinematic and photometric properties of these objects, in- 
cluding velocity dispersions, concentrations, stellar masses, 
light profiles, and phase space densities, are consistent with 
their eventual evolution into typical ~ elliptical galaxies. 

The coiTelat ions obeyed by th es e mergers an d rem- 
nants (e.g.. iGenzel et al.1 120011; iTacconi et alj 120021; 
iRothberg & Joseph! I2006albl: iDasvra et all 120071) are similar 
to e.g. the observed fundamental plane and Kormendy 
relations for relaxed ellipticals, and indicate that they will 
evolve onto these relations as their stellar populations age. 
This is further supported by the ubiquitous presence of fine 
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Structures such as shells, ripples, and tidal plumes in ordi 
nary ellipticals (e.g. Schweizer & Seitzer 1992; Schweize 
1996), whi ch are signatu r es of mergers (e.g . Ouinn il98 
.Hernquist & Ouinn! 119871: iHernquist & Spergel ,1992i) . and 
the clustering and mass density of ellipticals, which are 
consistent w ith passive evolution after formation in mergers 
dHopkins et a l. 2007c). These various line s of evidence are 
in accord with the " merger hypothesis" (iToomre & Toomr3 
Il972t Irbomrd [19771) . that elliptical galaxies originate from 
mergers of spirals. 

Numerical simulations performed during the past twenty 
years verify that major mergers of gas-rich disk galaxies can 
plausibly account for thes e phenom ena and elucidate the un- 
derlying physics. In Hopkins et alj (,2008d) . we describe the 
phases of evolution characterizing these mergers and their re- 
lationship to quasars, starbursts, and elliptical galaxies; we 
briefly summarize these phases here. Tidal torques excited 
during a m erger drive rap id inflows of gas into the centers 
of galaxies (iBames & Hernquist 1991, 1996), triggering star- 
bursts (Mihos & Hernquist 1996; Springel & Hernquist 2003j) 
and feeding rapid black hole growth tDi Matteo et aL,2005i) . 
Gas consumption by the starburst and dispersal of residual 
gas by s upernova-driv en winds and feedback from black hole 
growth dSpringel et al. 2005a) terminate star formation so that 
the remnant quickly evolves from a blue to a red galaxy. The 
stellar component of the progenitors provides the bulk of the 
material for producing the re mnant spheroid (Barnes 19881 
1 1 992!; iHernquislll 1 992l ll 99 3b') through violent relaxation, but 
the highest phase-space density material and central "cusp" 
come from dissipation in the final, merger-induced starburst. 

A major merger is generally required in order for the tidal 
forces to excite a sufficiently strong response to set up nu- 
clear inflows of gas and build realistic spheroids. Although 
simulations suggest that the precise definition of a major 
merger in this context is somewhat blurred by the degener- 
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acy between the mass ratio of the progenito rs and th e orbit 
of the interaction (lHernauisllll989l:lHernquist & Miho J 19951: 
Bournaud et alj I2005h . nu merical studies (lY ounger et aU 
2p07a b) and observations dPasvra et alJ l2006n Woods et ay 
2006) find that strong gas inflows and morphological trans- 
formation are typically observed only for mass ratios below 
^3:1, despite the greater frequency of higher mass-ratio 
mergers. In what follows, unless explicitly noted, we use the 
term "mergers" to refer specifically to major mergers (mass 
ratio < 3 : 1). 

Gas dissipation and star formation in the compact, merger- 
induced central starburst are key ingredients in these simu- 
lations and observations. It has long been recognized that 
many properties of ellipticals cannot be repro duced by purely 
dissipationless mergers of late-type galaxies ('Ostrikei) ll980l : 
ICarlberg 1986; Naab & Ostriker 2007 ): in particular, flieir 
high central phase-spac e densities and specific frequency 
of global clusters (e.g. ISchweizeJ Il998l) . degree of rota- 
tional support, a nisotropy, and minor axis rotation/kinematic 
misalignment (Barnes &Hernquisl[i99l iCox et alJ 120065 
|Naab et al. 2006), and isophotal shapes. However, these ar- 
guments do not apply when the merger constituents contain 
gas, which can radiate energy and hence increase phase-space 
densities (Lake 1989) and form new stellar populations. 

Observed absorption-line spectra and the red colors of el- 
liptical galaxies suggest that their stars formed at high red- 
shift (z > 1) and that little star formation has since occurred 
in them. If these galaxies were produced by mergers, the 
progenitors of many present-day ellipticals were high-redshift 
spirals. While relatively little is known about disk galaxies at 
high redshift, it is likely that they were more concentrated and 
more gas-rich than their low-red shift c ounterparts. Indeed, 
observational evidence (.Erb et'ani2006l) indicates that galax- 
ies at redshift z ~ 2 do have large gas fractions /gas '--^ 0.5, with 
some approaching /gas ^ 0.8-0.9. Thus, any attempt to un- 
derstand the formation, properties, and scaling relations of the 
present-day population of elliptical galaxies within the con- 
text of the merger hypothesis must consider gas-rich mergers. 

In pioneering studies employing simulati ons that incorpo- 
rated gas-dynamics and star formation, iMihos & Hernquis3 
(|l994a c) demonstrated that dissipation in disk-disk mergers 
leads to central starbursts similar to those observed in LIRGs 
and ULIRGs, leaving behind stellar populations that easily 
match the central phase space and surface densitie s of ob- 
served ellipticals. In particular, iMihos & Hernquis3 (Il994ah 
(hereafter, MH94) predicted that this process should imprint 
an observable signature in the surface brightness profiles of 
ellipticals, in the form of a steep departure from the outer 
Ide Vaucouleursid 19481) ^ r'/'^-law (which arises from violent 
relaxation of stars in the pre-merger disks) in the inner re- 
gions: i.e. the presence of a central cusp or extra light above 
the inwards extrapolation of the outer profile. Subsequent ob- 
servations (iHibbard & YunllT999t iRothberg & Joseohl 12004 
have now confirmed the existence of such excess light in at 
least a significant fraction of known gas-rich merger rem- 
nants. In principle, this could provide an indicator of how 
much mass was involved in the dissipational starburst, and 
therefore a critical constraint on t he formation history of the 
galaxy. iHibbard & Yuril (Il999l) : iRofliberg & JosephI (|2004 
also made first attempts to quantify the fraction of stellar 
mass or light involved; however, without a detailed compari- 
son with theoretical models to provide the appropriate "null" 
model, it is not clear how exactly to quantify the amount of 



"excess" in the light profiles. 

An estimate of how much gas is required t o match the cen- 
tral densities of ellipticals was provided by Hernqui st et al.l 
(11993), who used N-body simulations and analytic arguments 
to suggest that, without such a dissipational component, merg- 
ers would be deficient by ~ 15% of their mass relative to the 
inward extrapolation of an r'/"* law. iNaab & Truiillol (l2006h 
found a similar deficit by comparing a sample of surface den- 
sity profiles from collisionless (gas-free) disk merger rem- 
nants and observed merger remnants (but with reasonable 
agreement in profile shapes in the outer regions). Since true 
< ellipticals have a variety of profiles reflected in their Ser- 
sic indices, / oc exp{-(r/ro)'/"}, with n 2-4, this implies a 
range of extra light fractions ^ 1 - 20% needed to account for 
realistic populations of < L» ellipticals. 

At the time of the predictions made by MH94, these cen- 
tral cusps or extra light components, the relics of dissipational 
formation in a gas-rich merger, had not been observed. How- 
ever, with the advent of Hubble Space Telescope observations 
of the centers of elliptical galaxies, it is now well-established 
that typical < ellipticals exhibit such central light excesses, 
with mass ranges an d radii comparable to those originally 
predicted by MH94 ( Kormendv" 1 999'; Kormendy etalJlOoi; 
Ferrarese et al. 2006; Cote et al. 2007). The situation is some- 
what different for the most massive ellipticals, which appear 
to have central "cores" or "missing light," generall y accompa- 
nied b y slow rotation and boxy isophotal shapes dFaber et alj 
Il997h . It is believed, however, that this family of ellipti- 
cals has been further modified by subsequent re-mergers or 
"dry mergers" of two or more (initia l ly cus p y gas-rich merger 
remna nt) ellipti cals ("Guillen et al.l 120001; i Rest et all [20011: 
Ravin dranath et a l. 2001; Laine et al. 2003; L auer et al.ll2005l 
2007; Ferrarese et al.ii2006i:iCote et al. 2007; Kormendy et al.l 
2007). While the central profiles of these obj ects are certain ly 
interesting (and are considered in IHopkins et al.l (1200 8a')). 
they do not directly inform the history of gas-rich mergers 
and are therefore not the subject of our investigation here. 

In principle, there is great potential for understanding 
galaxy formation histories and structural properties by iden- 
tifying observed cusps with those predicted as the relics of 
gas-rich central starbursts by MH94. These cusps could en- 
able observers to determine individually the masses involved 
in violent relaxation versus dissipation in the merger, and 
to separate the stellar populations arising from both pro- 
cesses, further informing the origin of radial gradients in 
stellar age, metallicity, an d abundance in ellipticals (e.g. 
iMihos & Hernquisnll994bl) . The importance and fractional 
contribution of dissipation as a function of elliptical prop- 
erties such as mass and formation time could be assessed, 
and this could further be used as a direct test for the role 
of differential dissipation being a driver in e.g. the tilt of 
the fundamental plane, as suggested by the oretical models 
(iRobertson et al.1 l2006bt iDekel & Coxl l2006l) . Correlations 
with global kinematic properties can constrain and test the- 
ories, like those outlined above, that invoke dissipation to 
resolve long-standing objections to the merger hypothesis. 
The proposed links between ULIR Gs, quasars, and elliptical 
galaxy format ion in mergers (e.g. ISanders&Mirabellll996l: 
iHopkins et al.l l2006a b, 2008c, dj) similarly rely on dissipation 
occurring in mergers and, hence, anticipate that it should have 
a profound impact on the centers of ellipticals. 

Unfortunately, despite the major improvements in both ob- 
servations of the centers of elliptical galaxies and numerical 
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treatments of star formation and gas physics, there has been 
Httle attempt to use the extra Hght content of elHpticals as 
probes in the manner above. Rather, studies have mainly been 
limited to determining whether or not some extra light com- 
ponent is evident. This owes in large part to a lack of the- 
oretical underpinning: the original MH94 work simply pre- 
dicted that such cusps should exist; more detailed interpreta- 
tion and modeling was limited by spatial and numerical res- 
olution, and the inability to characterize gas physics with re- 
alistic empirically calibrated models for star formation and 
supernova feedback in the interstellar medium. There have 
been dramatic improvements in these areas in recent years 
(e.g. ' Springel & Hernquisj 120031: ISpringel et al.ll2005H) . but 

the issue has not been revisited. 

Here, and in a companion paper (i Hopkins et aT1 l2008b). we 
examine galaxy cusps or extra light components in both sim- 
ulations and observed galaxies, in order to identify the ori- 
gin of different components of elliptical surface density pro- 
files, their cosmological scalings and relevance for the forma- 
tion history of such galaxies, and their implications for global 
galaxy properties. In the present paper, we focus on the extra 
light in our simulations and in known gas-rich merger rem- 
nants, in order to isolate the relevant physics and to calibrate 
a means to recover the most important physical quantities in- 
volved. 

In §|2]we describe our suite of gas-rich merger simulations, 
and in § |3] we summarize the observational data set used in 
our comparisons. In §|4]we highlight the physical origin and 
nature of cusp or extra light components. In § |5] we com- 
pare different means by which to fit the surface density profile, 
and attempt to calibrate fitting methods in order to recover the 
physically distinct (dissipational versus dissipationless) com- 
ponents in merger remnants. In §|6]we briefly discuss the ef- 
fects of resolution and similar numerical issues involved, be- 
fore proceeding in §|7]to compare our simulations and fitted 
galaxy decomposition to observed systems. In §|8]we investi- 
gate how the amount of extra light depends on the gas content 
of the merging galaxies. In §|9]we use these comparisons to 
study how structural parameters of the outer stellar light and 
inner extra light component scale with e.g. galaxy mass and 
size. In §[TO]we calculate how these results and comparisons 
are affected by the choices of observed wavebands and time 
evolution after the merger Finally, in § [TT] we discuss our 
results and outline future explorations of these correlations. 

Throughout, we adopt a JIm = 0.3, ft\ = 0.7, i/o = 
70km s"' Mpc"' cosmology, and normalize all observations 
and models sh own to this cosmology. We also adopt a 
IChabried (120031) initial mass function (IMF), and convert all 
stellar masses and mass-to-light ratios to this choice. The ex- 
act choice of IMF systematically shifts the normalization of 
stellar masses herein, but does not substantially change our 
comparisons. All magnitudes are in the Vega system, unless 
otherwise specified. 

2. THE SIMULATIONS 

Our simulation s were perform ed with the parallel TreeSPH 
code GADGET-2 (Springel 2005), employing the fully conser- 
vative formulation ( S pringel & Hernquist 2002.) of smoothed 
particle hydrodynamics (SPH), which conserves energy and 
entropy simultaneously even when smoothing lengths evolve 
adaptively (see e.g., Hernquist 1993a; O'Sheaetal. 2005j). 
Our simulations account for radiat ive cooling, heating by a 
UV background (as in Katz et al.i 119961: iDave et al.1 1 19991) . 
and incorporate a sub-resolution model of a multiphase in- 



terstellar medium ( ISM) to describe star formation and su- 
pernova feedback (ISpringel & Hernquist! 1200 3h . Feedback 
from supernovae is captured in this sub-resolution model 
through an effective equation of state for star-forming gas, 
enabling us to st ably evolve disks with arbitrary gas fractions 
(see, e.g. Springel et al.ll2005bl: ISpringel & Hernquist! 120051: 
[Robertson et al.ll2006allcl) . This is described by the param- 
eter qsos, which ranges from ^eos = for an isothermal gas 
with effective temperature of 10'* K, to q'eos = 1 for our full 
multiphase model with an effective temperature 10^ K. We 
also compare with a subset of simulatio ns which adopt the 
star form ati on fee dback prescription from lMihos & Hernquist! 
( Il994a!lbll3. !1996 ), in which the ISM is treated as a single- 
phase isothermal medium and feedback energy is deposited 
as a kinetic impulse. 

Although they make little difference to the extra light com- 
ponent, supermassive black holes are usually included at the 
centers of both progenitor galaxies. These black holes are 
represented by "sink" particles that accrete gas at a rate M 
estimated from the local gas density and sound speed us- 
ing an Eddington-limited prescription based on Bondi-Hoyle- 
Lyttleton accretion theory. The bolometric luminosity of the 
black hole is taken to be Lboi = erMc^, where = 0.1 is the 
radiative efficiency. We assume that a small fraction (typi- 
cally « 5%) of Lboi couples dynamically to the surrounding 
gas, and that this feedback is injected into the gas as thermal 
energy, weighted by the SPH smoothing kern el. This fraction 
is a free parameter, which we determine as in Di Matteo et alj 
(2005) by matching the observed Mbh-o" relation. For now, 
we do not resolve the small-scale dynamics of the gas in the 
immediate vicinity of the black hole, but assume that the time- 
averaged accretion rate can be estimated from the gas proper- 
ties on the scale of our spatial resolution (roughly « 20 pc, in 
the best cases). 

The proge n itor g alaxy models are described in 
ISpringel et al.l (l2005bl) . and we review their properties 
here. For each simulation, we generate two stable, isolated 
disk galaxies, each with an extended dark matter halo with 
a Hernquist (1990 ) profile, mo t ivated by cosmolo gical 
simulations (N avarro et al J Il996t iBusha et af] l2005b . an 
exponential disk of gas and stars, and (optionally) a bulge. 
The galaxies have total masses Mvir = V^^JilQGHi)) for z = 0, 
with the baryonic disk having a mass fraction = 0.041, the 
bulge (when present) having nib = 0.0136, and the rest of the 
mass in dark matter The dark matter halos are assigned a 
concentration parameter scaled as in Robertson et al. (2006d) 
appropriately for th e galaxy mass and redshift following 
!BuIlock et al.! (!2001l) . We have also varied the concentration 
in a subset of simulations, and find it has little effect on our 
conclusions (because the central regions of the galaxy are, 
in any case, baryon-dominated). The disk scale-length is 
computed based on an assumed spin parameter A = 0.033, 
chosen to be near the mode in the A distribution measured in 
simulations dVitvitska et alj|2002l) . and the scale-length of the 
bulge is set to 0.2 times this. 

Typically, each galaxy initially consists of 168000 dark 
matter halo particles, 8000 bulge particles (when present), 
40000 gas and 40000 stellar disk particles, and one BH parti- 
cle. We vary the numerical resolution, with many simulations 
using twice, and a subset up to 128 times, as many particles. 
We choose the initial seed mass of the black hole either in 
accord with the observed Mbh-c relation or to be sufficiently 
small that its presence will not have an immediate dynamical 
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effect, but we have varied the seed mass to identify any sys- 
tematic dependencies. Given the particle numbers employed, 
the dark matter, gas, and star particles are all of roughly equal 
mass, and central cusps in the dark matter and bulge are rea- 
sonably well resolved. 

We consider a series of several hundred simulations of 
colliding gala xies, desc r ibed in detail in [Robertson et al.l 
(I2006blch and lCox etalJ (l20063bl) . We vary the numerical 
resolution, the orbit of the encounter (disk inclinations, peri- 
center separation), the masses and structural properties of the 
merging galaxies, presence or absence of bulges in the pro- 
genitor galaxies, initial gas fractions, halo concentrations, the 
parameters describing star formation and feedback from su- 
pernovae and black hole growth, and initial black hole masses. 

The progenitor galaxies have virial velocities Kir = 
55,80,113,160,226,320, and 500kms-i, and redshifts z = 
0,2, 3, and 6, and our simulations span a range in final 
spheroid mass Mbh ^ 10^- IO'^Mq, covering essentially the 
entire range of the observations we consider at all redshifts, 
and allowing us to identify any systematic dependencies in 
our models. We consider initial disk gas fractions (by mass) of 
/gas = 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, and 1 .0 for several choices 
of virial velocities, redshifts, and ISM equations of state. The 
results described in this paper are based primarily on simula- 
tions of equal-mass mergers; however, by examining a small 
set of simulations of unequal mass mergers, we find that the 
behavior does not change dramatically for mass ratios below 
about 3:1 or 4:1. This range is appropriate to the observations 
of merging galaxies used in this paper, which are restricted to 
major merger events. At higher mass ratios, our experiments 
([Younger et al. 2007a) suggest that gas can be channeled by a 
minor merger into the central regions of the galaxy, with an 
efficiency that declines in approximately linear fashion with 
the mass ratio. However, the resulting starburst forms a small 
bulge in a disk galaxy, so the decomposition into a Sersic pro- 
file representing the scattered stars and extra light component 
becomes ambiguous (both physically and observationally). 

Each simulation is evolved until the merger is complete and 
the remnants are fully relaxed, typically ~ 1 - 2 Gyr after the 
final merger and coales cence of the BHs. We then analyze 
the remnants following ICox et all (|2006b|) . in a manner de- 
signed to mirror the methods typically used by observers. For 
each remnant we project the stars onto a plane as if observed 
from a particular direction (we consider 100 viewing angles 
to each remnant, which uniformly sample the unit sphere). 
Given the projected stellar mass distribution, we calculate the 
iso-density contours and fit ellipses to each (fitting major and 
minor axis radii and hence ellipticity at each iso-density con- 
tour), moving concentrically from r = Q until the entire stellar 
mass has been enclosed. This is design ed to mimic obser- 
vational iso photal fitting algorithms (e.g. lBender et al.ll 19871 : 
lBender|[T98 8). The radial deviations of the iso-density con- 
tours from the fitted ellipses are expanded in a Fourier series 
in the standard fashion to determine the boxiness or diskiness 
of each contour (the 04 parameter). Throughout, we show pro- 
files and quote our results in terms of th e major axis radius. 
For further details, we refer to lCox et"an (1200 6b). 

We directly extract the effective radius as the projected 
stellar half-mass radius, and the velocity dispersion a as the 
average one-dimensional velocity dispersion within a circu- 
lar aperture of radius 7?e. Note that this differs from what is 
sometimes adopted in the literature, where /?e is determined 
from the best-fit Sersic profile, but because we are fitting Ser- 
sic profiles to the observed systems we usually quote both the 



true effective radius of the galaxy and effective radii of the 
fitted Sersic components. Throughout, the stellar mass 
refers to the total stellar mass of the galaxy, and the dynamical 
mass Mdyn refers to the traditional dynamical mass estimator 



Mdyn = k - 



(1) 



where we adopt k = 4.2 (roughly what is expected for a 
iHernquisl (|l990) profile, and the choice that most accurately 
matches the true enclosed stellar plus dark matter mass within 
in our simulations; although this choice is urelevant as 
long as we apply it uniformly to both observations and simu- 
lations). When we plot quantities such as R^, a, and Mdyn, 
we typically show just the median value for each simula- 
tion across all ^ 100 viewing directions. The sightline-to- 
sightline variation in these quantities is typically smaller than 
the simulation-to-simulation scatter, but we explicitly note 
where it is large. 

3. THE OBSERVATIONS 

We compare these simulations with t he sample of local 
remna nts of recent gas-rich mergers from lRothberg & JosephI 
( 120041) (henceforth RJ04). For these objects, RJ04 compile 
A'-band imaging, surface brightness, ellipticity, and a^/a pro- 
files, where the profiles typically range from ^ lOOpc to 
~ 10-20 kpc. These span a moderate range in luminosity (in- 
cluding objects from Mk ^ -20 to Mk ~ -27, but with most 
from Mk ^ -24 to Mk ^ -26) and a wide range in merger 
stage, from ULIRGs and (a few) unrelaxed systems to shell 
ellipticals. As demonstrated therein, these systems will gen- 
erally become (or already are, depending on the classifica- 
tion scheme used) typical ^ ellipticals, with appropriate 
phase space densities, surface brightness profiles, fundamen- 
tal plane relations, kinematics, and other properties. 

In terms of direct comparison with our simulations, the 
data often cover a dynamic range and have resolution com- 
parable to our simulations, and (experimenting with different 
smoothings and imposed dynamic range limits) we find it is 
unlikely that resolution or seeing differences will substantially 
bias our comparisons. We have converted the observations to 
our adopted cosmology, and compile other global parameters 
including e. g. kinematic proper t ies and luminosities in other 
bands from iRothberg & JosephI (l2006allbl) and the literature. 



We estimate stellar masses ourselves in a uniform manner for 
all the objects given the mean A'-band M /L as a function of 
luminosity from Bell et al. (2003), corrected for our adopted 
IMF. For objects where optical photometry is available, we 
have repeated our analysis using stellar masses derived from 
total K-hwA luminos ities and (^-V ) color-dependent mass- 
to-light ratios from Bell et al.' (2003) or from fitting integrated 
U BVRIJHK photometry of each object to a single stellar pop- 
ulation with the models of Bruzual & Charloti (120031) . and find 
this makes no difference to our conclusions. Of course, since 
these are young objects and we are not performing full SED 
fitting, there is some uncertainty in these estimates. We find 
a similar result, however, using d y namica l masses from the 
kinematics in Rothberg & JosephI (l2006al) - although these 
too are uncertain possibly owing to the incomplete relaxation 
of the systems. Within the considerable scatter in the observed 
properties, however, we find these details makes little differ- 
ence. 

Throughout, we will usually refer interchangeably to the 
observed surface brightness profiles in the given bands and the 
surface stellar mass density profile. Of course, stellar light is 
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not exactly the same as stellar mass, but in § [TO] we consider 
the differences between the stellar light and the stellar mass 
density profiles as a function of time, wavelength, and prop- 
erties of the merger remnant, and show that the /T-band results 
introduce little bias (i.e. are good tracers of the stellar mass); 
the Sersic indices and extra light fractions fitted to the A'-band 
profiles of the simulations are good proxies for the Sersic in- 
dex of the stellar mass profile and extra mass/starburst mass 
fraction, even close in time to the peak episode of star forma- 
tion. As has been noted in other works, many of these sys- 
tems have weak color gradients, empirically supporting the 
idea that there is generally only weak variation in M/L with 
radius. 

It is important to note that while we are not concerned about 
the absolute normalization of the profile (i.e. mean M/L), 
since we derive total stellar masses separately from the in- 
tegrated photometry, we must account for systematics that 
might be induced by change in M/L as a function of radius. 
The profiles in optical bands such as V require more care - 
when the system is very young (< 1 -2Gyr after the major 
merger-induced peak of star formation), there can be consid- 
erable bias or uncertainty owing to stellar population gradi- 
ents and dust. However, once the system is relaxed, the optical 
bands also become good proxies for the stellar mass distribu- 
tion. 

4. THE PHYSICAL ORIGIN OF "EXTRA LIGHT" 

We begin by considering a fiducial case which illustrates 
many of the physical processes relevant to the origin of ex- 
tra light in gas-rich galaxy mergers. Figure [T] shows the sur- 
face mass density of stars (from a random sightline) of a typ- 
ical, but extremely gas-rich, simulation. The merging galax- 
ies are similar to ~ disk galaxies (Sc-type galaxies) with 
a large initial gas fraction /gas ~ 0.8, and merge on an orbit 
representative of most random encounters. Despite the large 
gas content, the total stellar mass profile is quite smooth, at 
least down to the resolution limit of this particular simulation 
(^ 50 pc). However, there is a suggestion of a central extra 
light concentration or cusp within ~ 0.5- 1 kpc. (We choose a 
high gas fraction case merely to illustrate the salient features 
of our findings; as we detail in what follows, mergers with 
lower gas fractions yield qualitatively similar results.) 

The simulation allows us to decompose the stellar mass pro- 
file according to the time when various stellar populations 
were formed. We therefore consider three classes of stars. 
First, the stars that were present in the initial stellar disks at 
the beginning of the simulation ("disk stars") together with 
those stars formed before the final coalescence of the galax- 
ies and starburst phase ("pre-starburst" stars). Although the 
latter includes some stars that originate in a small starburst 
caused by the first passage of the two galaxies, most of these 
stars are formed quiescently in the two galactic disks before 
the disks are destroyed in the final merger. We identify these 
two simulated classes of stars as a single population because 
they are both "disk" stars entering the final merger, and it is 
of course arbitrary where (or at what time) we initialize our 
simulations, relative to the final merger, and therefore what 
fraction of the stellar mass forms "before" our simulation be- 
gins. As is clear in the figure, these stars are scattered dur- 
ing the merger and undergo violent relaxation, producing an 
r'/'*-like surface density profile. In detail, these stars relax 
to a nearly perfect Sersic profile with a best-fit Sersic index 
of tis ~ 4.0, which is shown in the right panel of Figure[T] As 
stars, they are coUisionless, and therefore remain at large radii 



and relatively low phase-space densities characteristic of the 
scale lengths and concentrations of the initial disks, with an 
effective radius R^. « 5 kpc. 

Second is the "starburst" component, defined somewhat 
loosely as those stars which form within 250 Myr of the 
final merger and coalescence of the two nuclei (taken for 
convenience as the merger of the two central black holes, 
but adopting another definition such as the peak in the 
star formation rate yields a similar result). This compo- 
nent is much more compact (Rg = 0.2 kpc), characteristic 
of the central starburst size scales in observ ed mergers and 
ULIRGs (e.g. lDo\ynes & Solomoni 1998; Br yant & Scovilld 
Il999t iGenzel et al.ll200U iTacconi et al...2002)) . Because our 
definition is based on only the time during which the stars 
formed, it includes some stars which originate in the merging 
disks (rather than the central, compact starbursts), which can 
be seen in the tails of this component of the surface density 
profile. However, these contribute negligibly to the surface 
brightness profile and total mass budget of the starburst com- 
ponent. The central component does not follow an r'/"* law, 
but is more similar to an exponential (n, = 1) profile. We will 
return to the details of the exact shape of this central light 
component in §|9l but in general expect this - the central star- 
burst is formed in situ by gas dissipation and star formation, 
and does not undergo violent relaxation to transform it from 
an exponential-like to an r'/"^ profile. 

Finally, because of the extremely gas-rich nature of this 
merger, there is some amount of cold gas remaining after the 
final black hole growth and feedback phase. This quickly 
settles into a rotationally supported disk, and forms an em- 
bedded disk of "post-starbu rst" stars (described in detail in 
ISpringel & Hernquistll2005h . Figure [T] shows that this can be 
well-fitted by a relatively small exponential disk (Rg ~ 3 kpc). 
Although this embedded disk is int eresting fr om the perspec- 
tive of the galaxy kinematics (see ICox et al.l '2006b). it con- 
tains only a tiny fraction ^ 4% of the stellar mass, and, be- 
cause this is concentrated near the total galaxy effective ra- 
dius, it adds negligibly to the surface brightness profile at all 
radii. 

We therefore confirm the findings of MH94, that gas chan- 
neled to the central regions of the galaxy forms stars in a com- 
pact starburst and leaves a central excess in the light profile 
above the extrapolation of the light of the old stellar popula- 
tion (scattered from the disks of the merging galaxies). How- 
ever, even in Figure [H where the extra light or starburst com- 
ponent constitutes a large fraction (-^ 30%) of the galaxy stel- 
lar mass, the overall light profile is quite smooth, and does 
not show an obvious sharp transition to the extra light com- 
ponent. This is unlike the earlier results of MH94, and we 
discuss some of the reasons for this in § |9] but they primar- 
ily relate to our improved resolution and treatment of the ISM 
relative to MH94, as well as more accurate time integration. 

5. FITTING TO RECOVER THE PHYSICALLY APPROPRIATE 
"EXTRA LIGHT" 

We would like to be able to use the surface brightness pro- 
file of the merger remnant to estimate the contribution of the 
extra light - ideally as a means of estimating, from the amount 
of extra light, how much of the galaxy mass was involved in a 
compact central starburst. However, as noted in § H] the light 
profiles in our simulation remnants are quite smooth, even 
where the extra light fraction is large. This makes identify- 
ing the extra light component a non-trivial procedure, which 
is sensitive to the assumptions made in fitting. 
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Fig. 1 . — Projected surface mass density of stars in ttie remnant of a highly gas-rich merger Left: Total profile (black) is shown with contribution from various 
physical components: stars in pre-simulation disks and those formed in the simulated disks before the final merger (disk/pre-starburst stars), stars formed in the 
compact, merger-induced starburst at final coalescence (starburst stars), and stars formed from gas which survives the merger (post-starburst stars). Right: Sersic 
profile fits to each of these components, labeled as "outer bulge," "cusp/extra light," and "embedded disk," respectively (sum shown in black). Mass fractions, 
Sersic indices («.,) and effective radii of each component are plotted (n, = 1 for the embedded disk and cusp/extra light component). 
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Fig. 2. — Upper Left: Surface mass density profile of remnant of a highly gas-rich merger (black), decomposed into stars formed prior to the final merger 
(which are then violently relaxed; red) and stars formed in the merger-driven dissipational starburst (blue). The Sersic fit to each component is shown, with the 
Sersic index of the pre-starburst component and mass fraction of the starburst component labeled. Upper Right: Two-component (Sersic plus cusp or extra light; 
where the cusp is fit with an exponential profile) fit to the total Hght profile, with the Sersic index of the outer component and mass fraction of the inner (n, = 1) 
component, and rms scatter ( A/i) about the fit. Lower Left: Single Sersic function fit to the profile. Lower Right: Cored Sersic function fit. The two-component 
fit as we parameterize it accurately recovers the Sersic profile of the violently relaxed component and mass fraction of the starburst component. The other fits 
give physically misleading results in this case. 
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We begin by reducing a galaxy surface density profile to the 
two most relevant physical components: the "pre-starburst" or 
"disk" stars, i.e. those formed in the rotationally supported 
disks before the final coalescence of the galaxies, and the 
"starburst" stars, formed in the compact, dissipational final 
starburst. We neglect the embedded stellar disks which can be 
formed by gas remaining after the merger; as demonstrated 
in Figure [T] these contribute negligibly to the surface mass 
density profile in even the most gas-rich merger remnants. 

Figure |2] shows the surface density of a fiducial gas-rich 
example, decomposed into a pre-starburst and starburst popu- 
lation. Motivated by observational methods and our results 
from § m we fit^ {upper left) a Sersic profile to the pre- 
starburst component, and find that it follows a nearly perfect 
r'/'* law (lis = 4.05). Loosely motivated by the shape of the 
starburst stellar profile at small radii and its origin in a gas-rich 
dissipational starburst, we fit the starburst stellar component 
to an exponential Sersic law (n, = 1). We emphasize that nei- 
ther the simulations nor observations herein resolve the inner- 
most regions of the central light where the choice of its Sersic 
index makes a significant difference, so we adopt this for sim- 
plicity as a standard choice in fitting the extra light component 
(we discuss this choice in greater detail below and §|9]l. This 
describes the starburst mass profile quite well where it is im- 
portant to the overall surface density. It does not capture the 
extended stars at > 3 kpc formed during this time, but as dis- 
cussed in § m we can safely ignore them. The exponential fit 
accurately recovers the total mass in the starburst component 
and its effective radius, and the shape of this profile where it 
is relevant to the total surface brightness profile. 

Rather than fit directly to the physical starburst and pre- 
starburst components, we now attempt to fit a two-component 
model to the observed quantity, the total surface brightness 
profile [upper right). We therefore fit the surface brightness 
profile to a double Sersic model, with an outer component 
reflecting the pre-starburst stars with a free Sersic index, and 
an inner component reflecting the starburst stars with a fixed 
Sersic index «j = L The total surface brightness profile is then 

/tot=/i/:,exp|-Z7i (^-^^| + 4exp|-/7„ (^^^ |, (2) 

where rsb and rb are the effective radii of the = 1 and 
free n, components (which we identify with the "starburst" 
and "old bulge" or "pre-starburst" components, respectively), 
Isb and lb are the corresponding normalizations, and n is the 
Sersic index of the outer bulge or pre-starburst component. 
The constant b„ is the appropriate function of n such that rsb 
and rb correspond to the projected half-mass radii. Note that 
changing our choice of = 1 for the inner Sersic component 
will systematically shift the mass balance between the inner 
and outer components, and lower the fitted outer Sersic index 
(typically increasing the inner mass by 30% and lowering 
the outer Sersic index An, ^ 0.25, respectively, if we adopt 
ris = 2). 

We have therefore experimented with a range of choices 
for the inner Sersic index: we find that fixing Wj = 1 recov- 
ers, on average (across our entire suite of simulations), the 

^ Formally, we consider our simulation profiles outside of some multiple 
~ 3 — 5 times the resolution limit, or with a seeing connection appropriate for 
the RJ04 sample (much larger than our resolution limits). We equally sample 
the profile in log r over a dynamic range extending to the largest radii in the 
RJ04 sample, and weight each point equally assuming an intrinsic ~ 0. 1 mag 
point-to-point variance in the SB profile (the typical magnitude of residuals 
fitting arbitrary splines to the data). We have varied these choices and find 
that our fits and conclusions are not sensitive to them. 



correct mass of the physical starburst component and Sersic 
index of the pre-starburst stellar component, whereas differ- 
ent choices for the inner n, give systematic offsets in these 
quantities. This is not to say that n, = 1 is most appropriate 
for every object/sightline or that it yields the most physical 
results in every case: however, if we are interested in recover- 
ing the physical extra light component in a statistical sense, it 
is most appropriate. Furthermore, even where we arrive at a 
similar result if we fit it with a free Sersic index, we find that 
this introduces considerably larger uncertainties when we fit 
observed systems (owing to the various parameter degenera- 
cies). 

Figure |2] shows the outcome of this fit. The resulting 
model of the surface density profile fits the simulation well, 
with a variance of only Ap ^ 0.09 mag arcsec"^ (assuming 
fi oc -2.5 log/tot). This is comparable to the point-to-point 
variance in the profile of this simulation if we fit an arbitrary 
spline to the data, and thus reflects a genuinely good fit. We 
have tested and found that the fit parameters are stable with 
respect to the dynamic range and error weighting, and that the 
residuals (typically of this magnitude, Ap^Q.l mag arcsec"^) 
are independent of radius (i.e. there is no systematic bias in 
the functional form). More important, this fit, despite hav- 
ing no direct information about the physical components into 
which we decompose the brightness profile, recovers almost 
exactly the appropriate parameters for both profiles. The best- 
fit Sersic index (n, = 4.04 compared to 4.05) and effective ra- 
dius of the "outer" or "bulge" component are almost perfect 
matches to those fitted directly to the pre-starburst stellar pop- 
ulation. Likewise, the inner or extra light component is a close 
match to the physical starburst component, and the fit recovers 
the extra light fraction accurately (32.4% compared to 32.5%; 
a much smaller difference than reasonable uncertainties in our 
physical definition of the starburst component). We have re- 
peated this fitting procedure for several hundred simulations, 
and find in almost all cases similar good agreement between 
the best-fit components and the physical decomposition of the 
galaxy stellar mass profile. 

However, there are a number of local minima in the fitting, 
and care must be taken to fit the most appropriate physical 
model. We have, in some sense, tuned our two-component 
parameterization (choosing the functional form and fixing the 
inner Wj = 1) so that it, on average, accurately recovers the 
distinct physical components in our simulated merger rem- 
nants. Other parameter choices may systematically fail (or 
fail in certain regimes of either high or low true starburst mass 
fraction) to recover a physically meaningful decomposition. 

For example, Figure|2]also shows the results of fitting a pure 
Sersic function to the entire surface density profile (lower 
left). There is a reasonable fit to the entire profile with a sin- 
gle Sersic index = 6.50, very different from the Sersic in- 
dex which describes either the pre-starburst or starburst light 
components. Likewise, if we consider the excess light to be 
that light in the real profile above the prediction of the best-fit 
Sersic model, we would infer only a tiny extra light fraction 
~ 0.5%, essentially consistent with no starburst component 
(given the fit uncertainties). Although the fit is technically 
worse, with variance Ap = 0.19, the difference is not dra- 
matic, and by many observational standards would be con- 
sidered a good fit. Clearly, however, the fitted results have no 
physical meaning in this case. 

The problem becomes even worse if we add a degree of 
freedom and fit a "cored Sersic" profile (lower right), of the 
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form 



/ = /' 



1 + 



j/a 



exp 



-bn 



(3) 



(e.g. ' Graham et al.ll2003h . where r/, is the core break radius 
within which the profile breaks to a power law of slope 7, 
re is the effective radius and n the Sersic index of the outer 
light profile, and a is a parameter describing how rapidly the 
break occurs. Figure |2] shows the results of this fit, which is 
again good in a purely statistical sense, albeit worse than our 
best fit Sersic+extra light fit (A/i = 0. 16 mag). However, here 
the fit parameters become even more unphysical - the best fit 
Sersic index is a very steep = 9.5 and one actually infers 
that the system is a cored galaxy, with missing light relative 
to the best-fit Sersic profile. 

The reason for these catastrophic failures is that the extra 
light component blends smoothly with the outer pre-starburst 
light profile. By increasing the central surface brightness, the 
extra light component makes the overall profile appear steeper 
(concave up in the /i-r'/** projection), owing to the rise at 
small r. However, the cusp itself does not continue to rise 
steeply inwards (in most cases), so after steepening the best- 
fit Sersic index to fit the outer part of the extra light compo- 
nent, one is often forced to infer the existence of a core in 
the central regions. Again, these fits have no direct physical 
meaning, but they are not terrible matches to the light profile. 
This emphasizes that whenever fitting a parameterized profile 
to the data, one must take care to adopt a physically motivated 
procedure. Fortunately, there are some indications, from the 
purely observational point of view, that the Sersic only and 
cored Sersic fits are inappropriate (when applied in this man- 
ner and into the central regions of gas-rich merger remnants). 

First, they begin to fail at large radii - however, this is where 
the true nature of the Sersic profile of the outer light compo- 
nent is most prominent, so any failure at large r should be 
especially worrisome. Furthermore, when we examine the 
kinematics (e.g. ellipticity, boxy/diskiness, rotation proper- 
ties) along the major axis, one can often see a transition oc- 
cur in these properties where the extra light begins to dom- 
inate the profile (see § |7]), whereas one would expect no 
such change if a continuous Sersic profile was the meaningful 
choice. Finally, when fitting a cored Sersic profile in an inap- 
propriate case, one often infers a large missing light fraction, 
relative to what is observed in genuine cor ed galaxies (e.g. 
massive, boxy, slow-rotating ellipticals; see iKormendy et alJ 
|2007). One typically expects missing light masses in real 
cores of < 1% the ga laxy mass, comparable to the central 
black hole masses (e.g. lMilosavljevic et alJl2002h . 

Note that we are not saying that a pure Sersic or cored Ser- 
sic profile is always an unphysical parameterization of the 
galaxy light distribution - however, for gas-rich merger rem- 
nants, we know in our simulations and have good reason to 
believe observationally that there is some excess light compo- 
nent. In these cases, the results of these fits are demonstrably 
unphysical and can be misleading. 

We caution that our adopted fitting functions and their 
physical interpretations are appropriate only for the galaxy 
spheroid, and should not be applied at radii where a galac- 
tic disk dominates the light profile. Implicitly, our choice of 
functional form also has the advantage of providing a check 
that stellar disks and distinct kinematic components do not 
contribute significantly to the surface brightness profile. Be- 
cause it includes a Sersic and an exponential component, if 
a traditional disk plus bulge decomposition of the light pro- 



file were required (if such a decomposition were the best fit 
to the observed profile), we would see this in the fit - our in- 
tended "extra light" component would in fact dominate the 
light profile at large radii, with a larger effective radius and 
substantially lower effective surface brightness than the fitted 
bulge. W e discuss these cases in more detail in lHopkins et alJ 
( l2008bl) . but note here that only a few of the simulated or 
observed systems show this behavior - in other words, they 
are primarily true ellipticals, not SO galaxies (which should 
not be surprising, since they are primarily remnants of ma- 
jor 1 : 1 mergers). In our simulations, we can confirm from 
the stellar kinematics that such fits are only preferred when 
indeed a rotationally supported disk is a significant fraction 
of the light at large radii - i.e. so long as we restrict our- 
selves to considering galaxy spheroids, there is no ambiguity 
that our fits are physically distinct from and statistically su- 
perior to simple disk plus bulge decompositions. Ideally, in 
the cases where a traditional disk/bulge decomposition is re- 
quired (i.e. there is a prominent large-scale disk), we would 
fit a three-component model: an outer exponential disk, and a 
two-component (dissipational and dissipationless) bulge. Un- 
fortunately, even with the ideal data quality in our simula- 
tions, the degeneracies involved in fitting for three indepen- 
dent components are large, and in cases where a disk domi- 
nates the profile at large radii, there is no robust "lever arm" on 
the outer bulge profile shape. We therefore generally exclude 
these cases from our comparisons (although we still can and 
do compare these profiles directly to observed systems), but 
emphasize that including or excluding them from our simula- 
tion sample or the observed RJ04 sample makes no difference 
to our conclusions. 

We also emphasize that although there are some superficial 
similarities between o ur adopted p arametric p rofile decompo- 
sition and that in e.g. ICdte et al] ([2006) and iFerrarese et al] 
(2006!), the two are in detail significantly different and ad- 
dress very different spatial scales and physical properties of 
the galaxies. Typically, the "outer profile" we refer to ex- 
tends to and beyond (in our simulations) the limits of our 
ground-based photometry, corresponding to physical radii of 

20- 100 kpc, and our "inner profile" refers to the resid- 
ual from a central starburst at scales where a significant frac- 
tion of the galaxy mass becomes self-gravitating (see § |9]l, at 
~ 0.5- Ikpc. We stress again that we are not resolving in- 
wards of the central ~ 30-50pc, and our modeling should 
not be extrapolated to within these rad ii without considerable 
care. In contrast, the "outer profile" in IFerrarese et al.l ( l2006h 
is based on the HST ACS profiles, which extent to outer radii 
^ 1 kpc, and their "inner profiles" typically dominate the light 
profile at very small radii 0.01 -0.02/?(, (-- 10-40pc for 
most of their sample). This is more akin to separating our 
"inner" component itself into multiple sub-components - i.e. 
a starburst stellar component that blends (as we have shown) 
relatively smoothly onto the outer, violently relaxed stars and 
an innermost nuclear component. The authors themselves ad- 
dress this, and denote these nuclear excesses as central stellar 
clusters. Such systems may indeed be present (and could be 
formed in the same dissipational starburst which we model): 
but if so they are distinct subsystems sitting on top of the star- 
burst light component, which we do not have the ability to 
model or resolve in our simulations or the observations (the 
typical nuclear star cluster scales are well below the resolu- 
tion limits of the observations from RJ04). Therefore, while 
the two approaches may yield complementary constraints and 
some similar conclusions, we caution that our results are not 
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Fig. 3. — Surface brightness profile as a function of spatial resolution in 
a resolution study. The profiles are converged outside of roughly a softening 
length. 

directly comparable and are specifically designed to trace dis- 
tinct physical structures. 

6. THE EFFECTS OF RESOLUTION 

Given the small physical sizes of the starburst or extra light 
components in the simulations, it is important to verify that 
finite spatial resolution is not biasing our conclusions. In gen- 
eral, we have experimented with a wide range of simulation 
spatial and mass resolution scales and particle numbers, and 
find that all of the results in this paper (except where other- 
wise explicitly noted) are robust to the effects of resolution. 
However, we briefly present a detailed resolution study to ex- 
amine these issues. 

Figure [3] plots the mean (averaged over ^100 sightlines) 
surface mass density profile in a series of otherwise identical 
simulations with increasing spatial resolution, where the grav- 
itational softening length decreases from w 150pc to 30pc. 
Note that the SPH gas smoothing length can be (and is) much 
smaller than this, especially in the dense, central regions of 
the galaxy during the starburst phase. 

Unsurprisingly, the profiles flatten within several softening 
lengths. However, they are well converged at larger radii (all 
the differences between the mean profiles plotted at > 5 soft- 
ening lengths are within typical sightline-to-sightline varia- 
tions). The primary difference seen in higher resolution sim- 
ulations is that the profile continues to rise to smaller and 
smaller radii. Note that there is still some well-resolved curva- 
ture to the profile - the central extra light or starburst compo- 
nent is not rising as steeply as a pure power-law, for example. 
However, we expect from this trend that the profiles, with fur- 
ther improved resolution, should continue rising to small radii 
< 30-50pc. Below these scales, we caution that, regardless 
of the numerical simulation resolution, our description of star 
formation and the ISM is no longer appropriate. Recall, we 
describe the ISM and star formation in a statistically volume- 
averaged sense, representative of a balance between the mul- 
tiphase components of the gas. By these smallest scales, we 
approach the physical sizes of molecular clouds and indi- 
vidual stellar clusters, and therefore cannot make physically 
meaningful claims without incorporating these gas phases and 
physics in a complete manner in the simulations. 

Given these issues, it is clear that simulation resolution does 
affect the shape of the central regions of the extra light com- 
ponents, and hence we cannot make strong conclusions re- 



0.35 
0.30 
0.25 
•I 0.20 
2 0.15 

LL 

0.10 

0.05 
0.00 

0.35 
0.30 
0.25 

.| 0.20 

o 

5 0.15 

LL 

0.10 

0.05 
0.00 



EL 



1 



3 



-0.1 



0.0 



0.1 

M, 



extra light 



0.2 

/M. 



0.3 



0.4 



Fig. 4. — Distribution (across ~ 100 sightlines) of outer Sersic index and 
inner/extra light mass fraction fitted to the resolution study profiles in Fig- 
ure|3](black, blue, and red for e, = 150, 70, 30 pc, respectively). The fits are 
reasonably well-converged even at relatively poor (150pc) spatial resolution. 

garding their innermost shapes. However, we also see that the 
outer profile (the pre-starburst component) is well-converged, 
and that the simulations exhibit a transition to an extra light 
component at similar radii for all resolutions. The resolution 
limits do not affect the radii where the excess light begins to 
depart from the outer profile (~ 0.5 - 1 kpc), or most of the 
range of the observed RJ04 profiles with which we compare 
(> lOOpc). Figure |4] shows the distribution of Sersic indices 
fitted to the outer component of these simulations (across all 
sightlines), and the distribution of the fitted extra light compo- 
nent fraction. Within the errors, these distributions are more 
or less fully converged even for our lowest-resolution simu- 
lation. Similarly, if we directly extract the "starburst" and 
"pre-starburst" components from the simulations, we find that 
their total mass fractions, effective radii, and surface bright- 
ness profiles (outside a few softening lengths) are converged. 
We also show in § |8] that the size-mass relation of the extra 
light components in the simulation is robust over the entire 
range of resolutions we have studied. We therefore conclude 
that resolution limitations, while constraining our ability to 
follow the profile of the extra light to the smallest radii, do 
not bias our overall estimates of the galactic structure or con- 
tribution of the extra light to the total mass. 

7. COMPARISON WITH OBSERVATIONS: THE EXTRA LIGHT IN 
LOCAL MERGER REMNANTS 

Bearing these cautions in mind, we apply our fitting proce- 
dure (calibrated from our simulations) to the sample of ^ 50 
local merger remnants from RJ04 which spans a representa- 
tive range in their post-merger stages (from LIRG and ULIRG 
systems still undergoing some star formation to shell ellipti- 
cals which are completely relaxed). 

We first illustrate the caveats from § |5] with a representa- 
tive example. Figure |5] shows the /T-band luminosity profile 
of UGC6, a merger remnant and LIRG in the sample of RJ04. 
RJ04 fit their profiles to single Sersic laws, and thek best fit 
is shown - they find a very large Sersic index «j = 10. RJ04 



10 



Hopkins et al. 




0.8 1.0 1.2 1.4 1.6 1.8 

Fig. 5. — Observed A'-band surface brightness profile (points) of the LIRG 
merger remnant UGC 6. Top: Single Sersic function fit (solid line) from 
RJ()4, with the variance about the fit A/x. Sersic index, and extra light frac- 
tion. Dotted vertical line marks the half-light radius. Bottom: Our two com- 
ponent fit (total, inner, and outer as solid, dashed, and dotted). The fit qual- 
ity is dramatically improved, and the outer Sersic index better describes the 
shape of the light profile at large radii. 

note that such a high Sersic value is probably unphysical, and 
parameters deriv ed from it should be r egarded with caution. 
Supporting this, iNaab & Truiillol ([2006) note that collision- 
less merger remnants, which can be reasonably approximated 
by a single Sersic law, do not result in such large Sersic in- 
dices. RJ04 also make a first attempt to quantify the extra 
light in their objects by comparing the total luminosity to the 
inward extrapolation of a r'/'^-law (different from the best- 
fit free Sersic index fit), similar in spirit to our approach but 
with less allowance for the detailed structure of the inner and 
outer portions of the galaxy. It is clear that the observed pro- 
file indeed shows considerable structure not captured in this 
fit (overall it is rather poor, with A/i = 0.24 mag), and that the 
profile at large r is in fact concave down in this projection (in- 
dicating Hs < 4), not concave up (Wj > 4). Furthermore, the 
status of this object as a LIRG indicates that the merger must 
have been highly gas-rich; such gas-rich events are expected 
to produce more disk-like outer light profiles, not large Sersic 
indices. 

Fitting instead to a Sersic-i-extra light component yields a 
much better fit in a statistical sense, even accounting for the 
added degrees of freedom (A/i = 0.03), as well as much more 
physical parameters - the outer Sersic index n, = 2.1 matches 
the concavity of the light profile and is expected in a gas-rich 
event, the profile is matched near the effective radius, and the 



shape of the extra light component is well-fitted. The result- 
ing extra-light component is much larger with this fit, ^45% 
(compared to ^ 16%) - this may seem high relative to what is 
typically quoted for such systems, but this has to do with our 
choice of fitting method. Indeed, such large fractional extra 
light components are reasonably common in simulations of 
very gas-rich mergers expected to produce LlRGs (see the ex- 
ample in § nil, and it is also observationally quite reasonable 
to see ~ half the ^T-band luminosity of a LIRG originate in 
a relatively compact ^ kpc nucleus (note though, as demon- 
strated in § [To] that the brightest LIRG and ULlRGs are the 
objects for which, even in K-bwd, there may be more extra 
light than extra mass, owing to stellar population effects). 

Figures l6lfT0l expand on these results by showing a com- 
parison with a representative subset of the RJ04 sample. (For 
comparison with the entire RJ04 sample, see §|A]and Table[T]) 
In each, we show the observed /T-band luminosity profile, 
with our best-fit two-component model, and note the fitted 
outer Sersic index, extra light fraction, and variance about the 
best-fit. We also directly fit our simulations to the observed 
profiles - i.e. considering the surface brightness profiles of our 
entire library of simulations, and find those which match the 
observed systems most closely. Because our simulations sam- 
ple a finite, discrete number of total galaxy masses, we allow 
the precise normalization of the simulation profiles to vary 
freely to match the observed data points (although we insist 
that the best-fit simulated galaxy mass be within a factor of 
~ a few of the estimated true galaxy mass, to make our com- 
parisons as physically robust as possible). We show the three 
best-fit simulations to each profile (with the variance of the 
data about the profile), along with the Sersic index directly 
fitted to the pre-starburst stellar component of each and the 
mass fraction in the starburst component. We plot these com- 
ponents for the best-fit case. Errors refer to the range of a 
given parameter (e.g. Wj of the pre-starburst component) in all 
simulations which are a comparably good match to the obser- 
vations. 

We also show 04 and ellipticity profiles of observations and 
simulations, although we emphasize that these are not fit- 
ted, and exhibit much greater sightline-to-sightline variation 
than the surface brightness profiles (thus should be considered 
somewhat less robust me asures of a gree ment). We caution 
that, a s demonstrated in ICox et all (12006 b) and iNaab et aLl 
( I2OO6I) . subtle features such as the isophotal shape are 
much more sensitive (compared to the surface brightness pro- 
file) to e.g. the orbital parameters and merger mass ratio. Fur- 
thermore, it is likely that seeing effects bias the observed el- 
lipticity and fl4 near the resolution limit, and this may be sig- 
nificant to radii as large as ^ Ikpc; we also note that if the 
central regions have higher light-to-mass ratios owing to a nu- 
clear starburst, they will appear more spherical than they truly 
are (although discrepancies in these shapes amount to small 
~ 0.1 magarcsec"^ effects in the brightness profile). For ex- 
ample, the observed ellipticity profiles here all go to zero at 
small radii <kpc, as do ground-based profiles of many nearby 
young or shell ellipticals - but in sufficiently high-resolution 
HST observations, it has been demonstrated that high ellip- 
ticities at ^ k pc scales continue (on average) at least down to 
< lOOpc (e.g. iLauer et al.ll2005l) . similar to what is predicted 
in our simulations. Nevertheless, as is also seen in lNaab et al.l 
(2006), the agreement between the observations of RJ04 and 
our simulations is quite good at most radii, even in the isopho- 
tal shape and ellipticity profiles. 

Both methods we have considered (fitting a parameterized 
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Fig. 6. — Observed, deconvolved major-axis A'-band light profiles of gas-rich merger remnants from RJ04 (points with error bars). Top: Fitted two-component 
model (inner, outer, and total as dashed, dotted, and solid lines), with outer Sersic index (n,), extra light fraction (/(,), and variance (Afi) shown. Middle: Direct 
comparison with simulated surface brightness profiles. The three simulations which match most closely are shown (red, green, and blue lines from best to 
third-best fit). Dashed line shows the physical starburst component of the best-fitting simulation. Range in outer «, and physical starburst mass fraction ( fsb) for 
the best-fitting simulations are shown, with the variance of the observations with respect to the simulated profile. Bottom: Observed disky/boxy-ness (a^) and 
ellipticity profiles, with the median (solid) and 25 — 75% range (shaded) con'esponding profile from the best-matching simulations above. Note that these are not 
fitted for in any sense. This panel shows the results for the ULIRGs, which require large gas contents and large extra light masses. 



profile and directly comparing to the simulations) have ad- 
vantages and disadvantages - a parameterized fit is of course 
independent of the simulations, and conclusions can be drawn 
from it even if the simulation physics are incomplete. On the 
other hand, there is always bias introduced in some parame- 
terized fit, and certain features in the light profiles (e.g. varia- 
tion at large radii from shells or ripples, and embedded disks 
at intermediate radii) can throw off our fitting routines, but 
are accurately reflected in the direct simulation mass profiles. 
Reassuringly, in almost all cases the parameters describing the 
best-fit simulations are the same (within the errors) as those 
from our parameterized fits. The variance about the best-fit 
simulations is also comparable to that from the fitted profiles - 
in a number of cases, the direct simulation profiles match the 
observed systems even better than a parameterized fit. This 
suggests that our fitting procedure is physically robust and that 
our simulations are reliably modeling real galaxy mergers and 
their extra light distributions. 

8. DEPENDENCE OF EXTRA LIGHT ON MERGER GAS CONTENT 

Having analyzed a large number of simulations and ob- 
served merger remnants, we can now study the properties of 
extra light or starburst components in a global sense. We ex- 
pect based on our above analysis that the extra light should 
be a reasonable proxy for the central, dissipational starburst 
component, and therefore should scale with the gas supply 
of the merging galaxies. Figure [12] shows this explicitly for 
a systematic survey of gas fractions. The mergers are other- 
wise identical, and similar to our ~ fiducial case in § ID 



but vary in initial gas fraction from /gas = 0.05- 1. We plot 
the sightline-averaged total surface brightness profile and me- 
dian best-fit extra light and outer/pre-starburst Sersic compo- 
nents in each case. It is clear that the extra light component 
systematically increases with increasing gas fraction, as ex- 
pected, typically containing ^ \/2 the initial gas mass (al- 
though this exact value depends on a number of other condi- 
tions, as shown below). The profile of the scattered, violently 
relaxed stars, on the other hand, is basically independent of 
gas fraction, as expected if these are simply being scattered 
from a rotationally supported disk to a partially violently re- 
laxed (Sersic-law n, ~ 2.5-3.5) object. 

Fi gures [T3] & [l4l extend this comparison to our entire suite 
of simulations and the fitted observations. We plot the ex- 
tra light fraction as a function of the stellar mass {upper left) 
or initial gas fraction of the simulation {upper right). As 
expected, the mean extra light fraction increases with initial 
gas content, but there is a significant residual dependence on 
the galaxy stellar mass and orbital parameters. At fixed ini- 
tial gas fraction, the extra light fraction declines with stel- 
lar mass, in a roug hly power-law fashion as /extra oc M^''-^^ . 
As demonstrated by [Robertson et all (l2006bl) . this is a conse- 
quence of the scaling of star formation efficiency with galaxy 
stellar mass. Higher mass galaxies have higher surface den- 
sities, and therefore exhaust their gas more rapidly - i.e. 
by the time the final merger and starburst occur, a typical 
^ IO'^Mq galaxy will have a smaller remaining gas frac- 
tion available to participate in the compact starburst than a 
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Fig. 7. — As Figure[6] Two LIRGs. which also require large gas content. 




Fig. 8. — As Figure[6] Two merger remnants with very prominent (deviant from the outer Sersic law) central excess light components. 
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Fig. 1 1. — As Figure |6] Two merger remnants from the study of lHibbard & YunI jl999h . The central gas identified in that work is shown (black solid line 
in middle). These are remnants, so the central excess/starburst component as we identify it (and in the con'esponding simulations) has already formed, and 
blends smoothly with the outer profile. The gas remaining corresponds to gas surviving the starburst (not part of the excess light), forming embedded disks (see 
Figure [T). 
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Fig. 12. — Median surface brightness profile (solid), with our inner (dashed) plus outer (dotted) component decomposition, of a typical simulation as a function 
of initial gas fraction. Outer Sersic index and extra light fraction for each are shown. All else being equal, the strength of the extra light component increases 
systematically with progenitor gas content. This does not necessarily mean that the excess light departs more shaiply from the outer profile shape. 
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Fig. 13. — Fitted extra light fraction/mass as a function of stellar mass {top left), initial (simulation) gas fraction (top right), gas fraction involved in the final 
coalescence/merger (bottom left), and total gas mass in the final merger (bottom right). Shown are the results for simulations: color encodes gas fraction (red, 
orange, green, cyan, blue, and black for initial /gas '^0.1, 0.2, 0.4, 0.6, 0.8, 1.0), symbol shape encodes orbital parameters (circles and triangles are typical 
random orbits, squares coplanar, stars polar), and symbol fill encodes initial redshift (filled for z < 1 typical progenitor disks/halos; open for ^ > 3 compact 
progenitor disks/halos). The plotted simulation values are the median across 100 sightlines; sightline-to-sightline differences are illustrated in Figure|4] 
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Fig. 14. — As Figure [T3l but with the simulations shown in black and the observed systems from the sample of RJ04 (magenta diamonds). Where multiple- 
component stellar population analysis allows an estimate of the mass fraction formed in the merger-induced starburst, we have plotted the observed systems in 
the lower panels. 
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typical M, ^ lO'^M© galaxy, even when both are initialized 
with the same gas fraction. Similarly, the orbital parameters 
determine how long it will take the galaxies to merge, and 
therefore how much of the gas will be consumed in the disks 
before the final merger event. 

The physically relevant quantity is therefore not the gas 
fraction at the beginning of our simulations (which, as noted 
before, is somewhat arbitrary), but the gas fraction of the sys- 
tem just before the final merger/coalescence of the galaxies. 
If we consider the gas fraction ^ 0.2 Gyr before the peak of 
the final starburst or coalescence of the black holes (roughly 
the "beginning" of the starburst, in most systems). Figure [T3] 
(lower panels) shows that there is a good, linear correla- 
tion between this pre-merger gas fraction and the extra light 
fraction in the remnant (/extra ~ /gas. pre-merger)- Accounting 
for this, there is no significant dependence on galaxy mass 
or orbital parameters. The only systems which deviate sig- 
nificantly tend to be pathological, either extremely gas rich 
/gas = 1, very high redshift z = 6, or e.g. perfectly co-planar 
mergers. In all of these cases, the systems can retain a large 
amount of gas eve n after the merger and fo rm relatively large 
new gaseous disks dRobertson et al.l2006ah . which are not ac- 
counted for in our fitting. Even in these situations, the devia- 
tions are not large - plotting the total mass in the extra light 
component against the gas mass present just before the final 
merger shows a tight correlation, as expected. 

Similarly, we find that changing properties of the simula- 
tions such as the presence or absence of an initial bulge, the 
concentration of the progenitor halos and disks, the presence 
or absence of a supermassive black hole, and the treatment of 
star formation and the ISM equation of state can all appear to 
change (albeit by a smaller amount than directly changing the 
gas fraction of the merger) the final mass in the starburst or 
extra light component. In all these cases, however, the effect 
is indirect - these choices influence how efficiently gas is con- 
sumed and/or expelled before the final merger, and therefore 
how much is available to participate in the starburst. For a 
fixed gas mass at the time of the final starburst, the extra light 
component mass is independent of these effects. 

Figure [14] summarizes the results of our fitting to the en- 
tire sample of RJ04. We note that one should regard the stel- 
lar mass estimates (based just on the /T-band luminosities) as 
somewhat uncertain - howe ver, using dynamical masses from 
iRothberg & JosepE (l2006ah yields a similar result and using a 
more sophisticated stellar mass estimator makes little differ- 
ence within the considerable scatter The observed systems 
occupy a similar locus in the /extra -Af* plane (upper left) to 
our simulations. There are a couple with very low /extra - 
however, without exception, these also have rather low best 
fit outer Sersic indices w 1 .0, and are probably better clas- 
sified as SO or spheroidal galaxies, for which our fitting rou- 
tines are not appropriate and the best-fit /extra is not necessar- 
ily physically representative. 

We cannot, of course, assign a meaningful initial gas 
fraction to the observed systems. However, it is in prin- 
ciple possible, by studying the stellar populations in suf- 
ficient detail, to estimate the mass fraction which formed 
in a recent, central starburst event (as opposed to the 
more extended quiescent star formation history). Unfortu- 
nately, there are still a number of degeneracies, and this re- 
quires deta iled obs ervations, but it has been a tt empted by 
iTitus et al.' (1997); Schweizer & Seitzer (1998); SchweizejJ 
71996); Schweiz er & Seitzer (2007); Reichardt et_al, (20011); 
iMichardI (l2006l) for several of the observed systems. Adopt- 
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Fig. 15. — Outer Sersic indices in observed and simulated gas-rich merger 
remnants, using our two-component decomposition. Points are as in Fig- 
ure [13] Gas-rich merger remnants have characteristic outer profile Sersic 
index n ~ 2 - 3, without a strong systematic dependence on mass or other 
properties. 

ing their estimates for the mass fraction in the newly formed 
stellar populations, we compare our fitted extra light fractions 
to the gas mass which participated in the starbursts in these 
objects. We find a similar tight correlation - although there 
are only a few objects for which sufficiently accurate stellar 
populations are available to allow this comparison, they all 
suggest that our fitted extra light component is indeed a good 
proxy for the mass fraction which was involved in the central, 
merg er-driven starb u rst. 

In 'HopkinsetiD (12008b'), we compare the properties of 
progenitor disks required to form realistic merger remnants 
and ellipticals to the observed properties of disks, and find 
good agreement. Here, we briefly note that the implied gas 
fractions in Figure [14] agree reasonably well with those ob- 
served in local disks of the same stellar mass (for stellar 
masses ^ 0.1 - IM*, where most of our sample lies, obser- 
vations imply typical disk gas fractions ~ 0.2 at z = to ^ 0.3 
at z = 1, with a factor ^ 2 intrinsic scatte r at both redshifts; 
iBeU & de Jonill2000t IShi)ley et al.ll2005h . Likewise, we ob- 
tain the fits to the observed surface brightness profiles in Fig- 
ures |6][TT]given progenitor disks that obey the baryonic Tully- 
Fisher and size-stellar mas s relations (observed to evolve only 
weakl y from z = 0--l; e.g . iBeU & de Jongi200U iFlores et alJ 
l2006tlTruiillo et al.ll20()6l) . We have performed some limited 
studies of initial conditions which do not resemble observed 
disks (e.g. initial disks that are very compact for their mass), 
and find that these fare poorly at simultaneously matching the 
observed light profile s and stellar masses, o r require unlikely 
gas fractions (see also lHopkins et al.]|2008al) . In short, the ob- 
served profiles are consistent with the expected remnants of 
typical local gas-rich mergers. 

9. STRUCTURE AND SIZE OF THE EXTRA LIGHT COMPONENT 

In Figure [15] we examine how the Sersic indices of 
the outer or pre-starburst light component vary with galaxy 
properties. Interestingly, there is no significant trend with 
galaxy mass. Thi s appea r s contrar y to th e conclusions of 
iPrug niel & SimienI (Il997h : iGrahamI (l2001l) : iFerrarese et aP 
(2006); however, we emphasize that we are fitting only gas- 
rich major merger remnants. Their samples include a number 
of spheroidal galaxies and "pseudobulges" at low masses (at 
the lowest values) and cored, boxy, slowly rotating ellip- 
ticals at the highest masses (and highest values). To the 
extent that gas-rich merger remnants dominate the population 
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at somewhat intermediate masses, the distribution of Sersic 
indices we reco ver is consistent with th eir estimated trends. 
Supporting this, iKormendv et alj (l2007b find a similar, rela- 
tively constant distribution of Sersic indices if they consider 
just the elliptical galaxies in their sample which have some 
central excess light. Furthermore, our fitting procedure should 
not be directly compared to that in lPrugniel & SimienI ( Il997h 
and others, because we fit an outer bulge and inner extra light 
component simultaneously, whereas they fit profiles to a sin- 
gle Sersic index. 

Physically, this independence of n, on other parameters is 
expected - the gravitational physics which scatter the stel- 
lar disks in violent relaxation are self-similar. It is only the 
gas physics of dissipation and star formation, responsible pri- 
marily for the extra light component, that break this self- 
similarity. This is explicitly clear if we c ompare our simu- 
lations to those in iNaab & Truiillol (l2006h . who find that in 
simulated collisionless (gas-free) disk merger remnants - i.e. 
systems for which the entire profile is by definition part of the 
"outer," violently relaxed component - there is no significant 
dependence of the Sersic index on mass, effective radius, or 
merger mass ratio. At the lowest masses, some of our simu- 
lated remnants have very low Sersic indices rig ^ I - these are 
generally low mass, extremely gas rich systems which form 
large disks after the merger, and should not therefore be con- 
sidered typical spheroids. 

We compare with the Sersic indices of the RJ04 sample. 
They occupy a similar locus and have a similar distribution 
to the simulations. There are a few objects with extremely 
low Sersic indices; as mentioned in §|8]and Naa b & Trujillol 
([2OO6), these are probably systems that are actually disk- 
dominated, and have small extra light components. For the 
systems that are genuinely bulge-dominated, the distribution 
of Sersic indices observed is consistent with the simulations. 
We also compare with the Sersic indices of relaxed ellipti- 
cals that have central cusps or extra light components, from 
iKormendy et al.l (l2007h . These trace an almost identical dis- 
tribution to the simulations and the RJ04 sample. Here, we do 
not see the anomalous objects with very low Hs because the 
systems are specifically selected to be bulge-dominated. In 
both the observed merger remnants and ellipticals, there is no 
trend of with galaxy ma ss, as predicted by the s imulations, 
a point discussed further in lHopkins et al.l (l2008bl) . 

Figure [16] shows the full distribution of Sersic indices in 
the simulation outer bulge or pre-starburst components. The 
distribution is smooth and can be approximated by a simple 
Gaussian with median Sersic index ~ 2.6 and Icr dispersion 
of An, = 0.75. The distribution of Sersic indices owes primar- 
ily to galaxy-to-galaxy variations: the figure shows the range 
of surface brightness profiles and fitted W j values for three ran- 
dom simulations, over ~ 100 sightlines. There is relatively 
little sightline-to-sightline variation in the surface brightness 
profile, and in general a small dispersion in Sersic index (dis- 
persion of Aris ~ 0.1-0.3 across viewing angles). Again, this 
is similar to the simulated di stribution in collisionless disk- 
disk merger remnants in Naa b & Truii llo (2006), which in- 
clude a wider range in e.g. merger mass ratios, but do not 
include gas and therefore should be strictly self-similar. 

For comparison with what is typically measured. Figure [TT] 
shows the distribution of Sersic indices measured if we simply 
fit a single Sersic profile to the entire surface brightness pro- 
file of each of the simulated and observed merger remnants. 
Fitting a single Sersic profile, we see that the distribution of 
Sersic indices is systematically pushed to higher values, with 
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Fig. 16. — Top: Surface brightness profiles of three random (~ L*, /gas ~ 
0.2 — 0.4) simulations with slightly different n, values. Solid lines show the 
median profile over ~ 100 sightlines, shaded range shows the 25—75% range 
across sightlines. Middle: Distribution of fitted outer h, values for each of 
the simulations. Bottom: Distribution across all sightlines and simulations of 
outer Sersic index in our simulated gas-rich merger remnants. A Gaussian fit 
to this distribution is also shown. 



a median n, ^ 4.0-4.5 for both the simulations and obser- 
vations, and a much larger "tail" towards high Sersic indices 
rts > 5. Fitted in this way, Sersic index does appear to depend 
(albeit weakly) on galaxy mass, as well as effective radius 
and surface brightness, similar to what has been seen in sam- 
ples of ellipticals (when fit in this manner) as discussed above. 
We caution, however, that as demonstrated in Figure |2l these 
fits are not physically representative of the outer light profile 
or central extra light. The presence of an extra light compo- 
nent makes the total central light profile steeper, leading to the 
higher ^ 4 typical of ellipticals that are fit to single Sersic 
indices - as compared to the ~ 2 - 3 that genuinely char- 
acterizes the outer, violently relaxed stellar populations in our 
simulations. Because, as we have shown, this outer profile is 
roughly self-similar, the dependence of the single Sersic index 
on mass and other galaxy properties seen here owes to differ- 
ences in the typical masses, shapes, and spatial extent of the 
extra light components, not to differences in the outer profiles. 
Even allowing for a fit to the full profile, RJ04 found (similar 
to what we see here) that the correlation between luminosity 
and Sersic index had a large scatter, and the correlations be- 
tween e.g. size and Sersic index were even less significant. 
In short, the behavior of a single Sersic law fit in the simula- 
tions and observations is similar to what has generally been 
observed in elliptical galaxies, but when fit in this manner this 
"typical" behavior reflects a combination of the inner extra 
light structure (which depends on the degree of dissipation 
and therefore other galaxy properties) and outer, violently re- 
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Fig. 17. — Top: Best-fit Sersic indices ob taine d from simulated and ob- 
served RJ04 merger remnants (style as Figure [TH when fit to a single Sersic 
profile (as opposed to our two-component fit). Bottom: Distribution of these 
single Sersic profile indices (histograms with and without error bars are for 
the RJ04 observed sample and simulations, respectively). Fitting to a sin- 
gle Sersic law yields higher Sersic indices and shows a dependence of ji, on 
galaxy stellar mass, radius, and surface brightness similar to what is "typ- 
ically" observed. As shown in Figure |2] however, this does not reflect the 
true shape at large radii (and should not be used to estimate e.g. the true ex- 
tra light fraction), but rather indirectly reflects the contributions of extra light 
changing the profile shape at small radii. 

laxed stars (which are approximately self-similar), rather than 
any single, robust physical component of the galaxy. 

We can, in principle, also attempt to measure the Sersic in- 
dices specifically of the cusp or extra light components. Un- 
surprisingly, they typically show disky n, 1 - 2 profiles. 
However, we caution that the inner regions of the cusp are 
strongly affected by our resolution limits. It is guaranteed 
that near the resolution limit the surface density profile will 
become flat, which results in a lower Sersic index. Therefore, 
although the general trend that the extra light component tends 
to have a diskier profile than the outer light appears to be true 
in the observations as well, the structure in the inner regions 
of the simulated extra light components should not be consid- 
ered robust. As discussed in § |5] our choice to fit the inner 
extra light component to an Wj = 1 law is not motivated by the 
actual shape at small radii (well inside the effective radius of 
the extra light), which neither simulations nor observations re- 
solve (and high-resolution observations of ellipticals suggest 
can have irregular structure owing to e. g. star clusters and the 
effects of black hole dynamics; see e.g. lLauer et al.ll2005l) : in- 
stead, we chose this because it fits reasonably well to where 
the extra light maps on to the outer profile and provides the 
most reliable recovery of the true physical starburst compo- 
nent. 

More robustly, we can study the effective radii of the extra 
light or starburst components in the simulations and observa- 
tions. Figure [T8]compares the spatial extent of the extra light 



component to the properties of the galaxy. The extra light 
component appears to be fractionally smaller in more mas- 
sive galaxies (relative to the effective radius of the galaxy), in 
both simulations and observations. Looking at this in detail, 
it is clear that this is a consequence of the trend in Figure [T3] 
- more massive galaxies tend to have less mass in their ex- 
tra light components. The extra light itself appears to follow 
a rough size-mass relation, shown in Figure [18] of the form 
Re, extra light oc M^^^^ Notc that by size of the extra light 
component, we refer explicitly to its half-mass projected ef- 
fective radius: this is not necessarily the same size scale at 
which the extra light component begins to dominate the total 
surface density profile of the galaxy. The appropriate compar- 
ison is with the effective radii of the extra light components 
fitted to the observations - where we see good agreement. By 
mass of the extra light component, we simply refer to our fit- 
ted extra light fraction times the total estimated stellar mass 
of the galaxies. 

Considering the evolution of the starburst component in the 
simulations, it appea rs to follow a general for mation mecha- 
nism, as outlined in Mihos & Hernquisj (11996 ). In the final 
stages of the merger, the gravitational torques acting on the 
gas remove its angular momentum, and it essentially begins a 
free-fall towards the central regions of the galaxy. The cool- 
ing time of the gas is almost always shorter than the dynamical 
time at this stage, so the process is not pressure supported and 
is primarily a gravitational collapse. This continues until the 
collapsing gas becomes self-gravitating: at this point, the gas 
shocks and establishes a quasi-hydrostatic central mass con- 
centration, and becomes stable against further collapse. The 
stability criterion is simply that the effective equation of state 
of the gas be 7 > 4/3, which is true at this time because the 
gravitational compression of the g as is effectively adia batic. 
Furthermore, in the classical McK ee & Ostrikei] (11977 ) pic- 
ture of the multiphase ISM, the densities are sufficiently high 
that the system is subject to thermal instability. In this state, 
the hot and warm diffuse phases of the ISM can radiate and 
cool efficiently into cold, dense clouds, which then form stars, 
whose supernovae and stellar wind feedback reheat the dif- 
fuse ISM. The important point is that, in this regime, cooling 
does not allow the system to lose effective pressure support, 
because the cooling of the diffuse ISM is balanced by energy 
injection from star formation and supernovae in the cold phase 
of the ISM. Regardless of the cooling time, then, one expects 
that the system cannot collapse further once it becomes self- 
gravitating, and will stall at this radius and continue rapidly 
forming stars until the gas supply is exhausted. 

Studies of the central regions of starburst systems sup- 
port this picture - typically, observations find a dense con- 
centration of molecular clouds forming stars at a high rate, 
with very high effective temperatures and pressures of the 
cumulative (diffuse) ISM gas (e.g. [ Sanders & Mirabelll996t 
ISolomon et all [19971 iBrvant & Scoville..l999f r Even if the 
system can contract, it does so on the local dynamical time, 
which, for a self-gravitating system, is comparable to the star 
formation timescale if the system obeys a Kennicutt-Schmidt 
type star formation law (iKennicuttil 1 998) . and so the system 
will exhaust gas in star formation in a self-similar manner as 
it contracts. 

The expectation of these physical models is then that the ef- 
fective radii of the extra light components will be determined 
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Fig. 18. — Top Left: Ratio of the half-light radius of the fitted e xtra light component to that of the fitted outer Sersic component, as a function of galaxy 
stellar mass. Points show simulations and observations, as in Figure [T3] Top Right: Size-mass relation of the extra light. Bottom Left: Potential of the extra 
light component (here defined as GM/R) as a function of that of the entire galaxy. Dotted line denotes a linear relation (indicating that the extra light is 
self-gravitating). In this space, there is no systematic dependence of extra light comp onent siz e on e.g. galaxy mass, orbital parameters, or merger redshift. 
Colors here denote the effective ISM pressurization and equation of state: red for the full lSpringel & Hernquist ( 2003) model, blue for a nearly isothermal model. 
Bottom Right: Same, but colors denote simulation spatial resolution, from ~ 300 pc to ~ 30 pc. Stars denote the resolution study simulations from Figure[3] 



by the condition that they become self-gravitating: 
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where a 1 is a constant which depends on the exact shape of 
the galaxy mass profile and stellar-to-dark matter mass ratio 
as a function of radius. Figure [18] plots this relationship, for 
both the simulations and the observed sample of RJ04. We 
find that indeed, the objects follow a tight, linear correlation 
as predicted, with a best-fit normalization awl, independent 
of the absolute value of galaxy mass, orbital parameters, or 
the initial redshift of the merging galaxies. 

If the spatial extent of the extra light component were not 
set by the self-gravitation condition, but by, for example, the 
pressure support of the ISM or some competition between the 
star formation and cooling timescale, one would expect the 
size to be systematically sensitive to the treatment of star for- 
mation and supernova feedback and the corresponding ISM 
equation of state. Figure [18] compares the size-mass relation 
for two values of our equation of state parameter ^eo,s, de- 
scribed in § [2] Our fully multiphase model for the ISM, with 
a high effective pressure, corresponds to <7eo,s = 1, and a model 
much closer to a completely isothermal ISM to ^eos = 0.25. 
Both cases follow an identical size-mass relation, suggest- 
ing that it is indeed simply the self-gravitation condition that 
sets the size of the starburst component. Note that the abso- 
lute sizes of the extra light or starburst components tend to 
be smaller in the ^eos = 0.25 case. This is because the more 
isothermal equation of state allows star formation to proceed 



more efficiently in the pre-merger disks, exhausting more of 
the gas supply before the final merger and leaving an extra 
light component of lower total mass. But, for fixed gas mass 
just before the final merger stages or (equivalently) fixed star- 
burst mass fraction, the two equations of state yield extra light 
components of the same physical size. 

It is also natural to wonder whether or not the sizes of the 
cusps are sensitive to our simulation resolution. Figure [18] 
compares the points on this effective size-mass relation as a 
function of the simulation resolution. From our lowest res- 
olution simulations with gravitational softening length « 
300 pc, to our highest resolution simulations with « 30 pc, 
the systems obey the same correlation. This is reassuring in 
that our resolution limits, while such that we do not resolve 
the inner regions of the extra light, do not cause us to system- 
atically overestimate the size of the extra light distribution. 
We have also tested a series of simulations of varying integra- 
tion accuracy, and find similar results. 

10. STELLAR POPULATIONS AND EVOLUTION OF EXTRA LIGHT 
COMPONENTS WITH WAVELENGTH AND TIME 

Figure [T9]plots the evolution of the surface brightness pro- 
file of a highly gas-rich merger (initial /gas = 0.8) remnant 
with time, in the observed B-band and /T-band, ignoring (for 
now) the effects of dust obscuration. We calculate the ob- 
served surface brightness profile in t hese bands using the stel - 
lar population synthesis models of Bruz ual & Char iot (2003"), 
given the ages and metallicities determined self-consistently 
for the stellar particles formed in our simulations, assum- 
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Fig. 19. — Surface brightness profile of an extremely gas-rich (/gas = 0.8) 
merger remnant in B-band (left) and A'-band {right), as a function of time 
after the coalescence of the two merging nuclei (colors, as labeled). Top: 
Absolute values, which show fading owing to passive evolution. Bottom: 
Profiles rescaled to the same value at ~ 30kpc, to show the change in profile 
shape. There is Httle evolution in A'-band, but substantial time dependence 
in B-band. In this case the system forms an embedded disk (around t 
0.5 Gyr) from gas which survives the merger, which is prominent in B-band 
(the "bump" at ~ 4 — 15 kpc) owing to its young age but negligible in A'- 
band (it contains only ~ 4% of the galaxy stellar mass). The young stellar 
populations in the embedded disk have faded in B-band by ^ 3 Gyr. 

ing a IChabrieil (l2003h IMF (although this choice primarily 
affects only the total luminosity, not the shape of the light 
profile). The initial stars in the simulation disks (in this 
case a relatively small fraction 1 - /gas = 0.2) are given pre- 
merger age, metallicity, and enrichment distributions appro- 
priate for the disk stellar mass according to the median best- 
fit T model star formation histories fit to local observed disks 
in iBell & de Jond (|2000') (thereby are guaranteed to lie on 
the observed star-forming galaxy mass-metallicity relation). 
We discuss the issues of stellar popu lation gradients in much 
greater detai in Hopkins et alJ (l2008 b). but note that we have 
also experimented with other stellar population models or fit- 
ted star formation histories, and with including initial gra- 
dients in the disk stellar populations, and find these choices 
make little difference (the central age, metallicity, and gradi- 
ents of interest are largely set by the dissipational starburst). 

The figure demonstrates that, although there can be signifi- 
cant overall fading in the near IR with time after the merger (in 
this case ^ 2 magnitudes over ~ 3 Gyr), especially for such a 
gas rich system where a large fraction of stars are formed in 
the final starburst, this has little differential effect with radius. 
The shape of the surface brightness profile in ^T-band is nearly 
time-independent, after the final coalescence of the two galax- 
ies. In optical bands such as B, however, there can be a much 
more dramatic time dependence. In this case, not only does 
the nucleus fade considerably relative to the outer profile, but 
one can see a prominent embedded disk appear at ~ 0.5 Gyr 
and fade out by ~ 3 Gyr after the merger. 

Because the system is so gas-rich, a significant quantity of 
gas survives the merger and cools to re-form an embedded 
disk. This process takes ^ a couple x 10** yr after the coales- 
cence, hence the disk is not as prominent immediately follow- 
ing the merger By mass, the disk represents only ^ a few per- 
cent of the galaxy stellar mass, similar to the case in Figure[T] 
However, because its stellar populations are very young, it 
appears prominently in B-band, constituting ~ 30-40% of 
the total light and making the system appear to be a clear SO 



rather than a true elliptical. By ^ 3 Gyr after the merger, the 
embedded disk (while still forming stars at a low rate) has 
aged sufficiently so that its mass-to-light ratio is not much 
different from the rest of the galaxy, and (given its small mass 
fraction, which is concentrated near R,,, where most of the 
spheroid mass also lies) it disappears again. Note that, to il- 
lustrate this, we have ignored dust, which will preferentially 
lie in the gas-rich star forming region of the disk, and make 
the difference between B and K bands less obvious. Still, the 
presence of such disks in optical (and not NIR light) is not 
uncommon in gas-rich merger remnants: the profile and evo- 
lution above is similar to th at seen (and inferred) for NGC 34 
(ISchweizer & Seitzedl2007h . 

The example above and observed systems emphasize that, 
at times close to the merger, the /T-band is a much more robust 
tracer of the stellar mass distribution. Our experiments sug- 
gest that we can have reasonable confidence that the /iT-band 
profiles of the observed systems will not change substantially 
with time - therefore our results would likely be unchanged 
if we analyzed the merger remnants at a time when they were 
more relaxed. Of course, then, it would not be clear that the 
systems are indeed gas rich merger remnants. The most sub- 
stantial change would probably be a reduction in the scatter 
about the mean trends defined in Figures [T3] and [TSl we typi- 
cally see little evolution in the median surface brightness pro- 
files after the merger, but unrelaxed features can cause con- 
siderably larger sightline-to-sightline variation. 

Figure |20] shows the evolution with time of the best- fit sur- 
face brightness profile parameters for two more representative 
cases, one with a large initial gas fraction and one with a small 
initial gas fraction. Unlike Figure [19] which was just a quali- 
tative illustration, we include dust obscuration by calculating 
the column dens ities through the gas in our simulations in the 
manner of Hop kins et alJ (l2005allbl) . We show the evolution of 
the absolute effective radius (i.e. independent of any fitting), 
the best-fit outer Sersic index, and the best-fit extra light frac- 
tion with time after the final coalescence, in terms of stellar 
mass, /T-band light, B-band light, and B-band light in the ab- 
sence of dust obscuration. 

The effective radius and extra light fraction in terms of stel- 
lar mass show little evolution with time. In other words, the 
profile is relaxed out to ^ in a short time (~ the dynam- 
ical time). The behavior in Wj, however, is more sensitive to 
the outer regions where the dynamical time is long, and there- 
fore more varied. We see in Figure |20] that, in one case, 
declines slowly with time - this is caused by a large sloshing 
from the merger, which moves stars to large radii and forces a 
fit with a more extended envelope (higher n^), until the slosh- 
ing relaxes. In the other case, there is a clear ringing or os- 
cillation in ris- This is caused by material on nearly radial 
orbits creating shells, as d escri bed in detail in Ouinn ( 1984); 
iHernqu ist & Ouinn (1987) and iHernquist &"Sperge]| 019921) . 
The oscillations occur as the material in shells moves in and 
out to large radii. As more orbits are completed, the shell 
stars phase-wrap, smoothing their projected density distribu- 
tion, and so the ringing in n, damps out. These behaviors are 
not unusual, but there is considerable variation in the effects 
on n, with time - it is difficult to define a typical behavior for 
the relaxation of n,. It is also clear that the degree to which 
might be different after the system is relaxed to that measured 
in the recent merger remnant stage is highly variable - there 
is no well-defined systematic bias (or even a typical direction 
of such a bias, to higher or lower n,) in at early times ver- 
sus late times. More likely, we expect that there will be less 
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Fig. 20. — Evolution of surface brightness profile parameters as a function of time after the coalescence of merging nuclei, for a highly gas-rich major merger 
(top) and a less gas-rich major merger (bottom), in terms of the stellar mass profile, observed X^-band (including stellar population age and metallicity effects 
and dust obscuration), observed B-band, and B-band ignoring dust obscuration (as labeled). Values shown at each time are medians across ~ 100 sightlines, the 
sightline-to-sightline variance is generally slightly smaller than the offsets in the plotted curves. Left: Half-light radius (directly from the light profiles). Middle: 
Outer profile Sersic index . Right: Extra light fraction. The evolution in Jij is largely real (and does not generally follow systematic patterns), and owes to 
relaxation of shells and other features; other evolution owes to stellar populations and dust. 



scatter in ris in the merger remn ant sample of RJ04 ori ce the 
objects have all relaxed (see also' Naab & Truiinoll2006l) . 

The /T-band profile reflects these trends - there is some 
(generally weak, even in the extreme high /gas case) time evo- 
lution owing to the weak dependence of K-hand M/L on stel- 
lar age. The sense of this is what is generally expected - since 
stars in the center of the remnant are formed in the final star- 
burst and younger, they have larger L/M ai early times, and 
more of the light appears to come from inner radii, biasing 
Re to slightly smaller values. The effect is generally appli- 
cable for only ~ a couple x 10^ yr after the merger, however, 
considerably shorter than the likely ages of most of the ob- 
jects in the observed sample (excepting, perhaps, those still in 
the ULIRG/LIRG class). After this time, the age difference 
is small relative to the absolute ages of the populations, and 
makes little difference to M/L as a function of radius. 

That is not to say, however, that the /T-band is free of bias. 
The measured quantities reflect the stellar mass profile with 
little additional time dependence, but are systematically off- 
set by a small amount. This is because both of the examples 
shown have strong metallicity gradients, with higher metal- 
licities in their central regions. In /T-band, L/M rises with 
metallicity, so this bias goes in the same direction as the age 
bias, but is more permanent. The metallicity gradient does 
not fade in the same manner as an initial age gradient. There 
is therefore a systematic bias towards underestimating Re and 
overestimating /exa-aiight by roughly ^ 20%. The magnitude of 
this bias, however, is directly sensitive to the metallicity gra- 
dient, and observed ellipticals tend to have gradients ranging 
from nil to comparable strength to those above (discussed in 
detail in iHopkins et al.l ( l2008bl) ). so the bias will range from 
object to object in this range. It is therefore difficult to sys- 
tematically correct our analysis of the observed systems - but 
even if we uniformly apply the most extreme correction for 
this fact, it makes no difference to our conclusions, and the 
offsets implied are considerably smaller than the object-to- 



object scatter in these parameters. In fact, considering an en- 
semble of simulations, one is far more likely to introduce bias 
in Re and /extra light by not including sufficient dynamic range 
in the surface brightness profile used for fitting. 

In B-band, the comparison is roughly as expected. If one 
ignores dust, there is an extremely strong time dependence of 
all the fitted parameters until ~ 1 Gyr after the merger. Dust, 
however, negates much of this, because the youngest stellar 
populations introducing the greatest bias also tend to be the 
most obscured. To the extent that this is typical, the biggest 
effects in B-band are also in the first ^ couple x 10^ yr, where 
stellar populations are young and substantial dust exists which 
can make the best-fit parameters highly variable and sensitive 
to viewing angle. Interestingly, at late times after the merger, 
the optical light profile reflects the stellar mass profile with 
less bias than the NIR profile. This is because, in B and other 
optical bands, the trend of L/M with metallicity is opposite 
that with age: L/M decreases with metallicity. Therefore the 
effects of age gradients and metallicity gradients actually tend 
to cancel one another out at late times. This suggests that, at 
early times after a merger (typical of the systems in the RJ04 
sample studied here) the NIR is a better tracer of the surface 
brightness profile, considerably less sensitive to age effects 
and dust-induced sightline-to-sightline variation. Once the 
systems are fully relaxed, however (provided there has been 
no new infall of gas and dust or new star formation event), op- 
tical light becomes similarly robust, and can actually trace the 
mass distribution with slightly less systematic bias. In prac- 
tice, however, even evolved ellipticals often have dust lanes 
or central dust concentrations, so care must be taken regard- 
ing their effects. 

11. DISCUSSION 

We have studied the origin and properties of extra light 
or cusps in mergers of gas-rich galaxies using a large 
suite of both numerical simulations and local observed 
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merger remnants. We con firm the original prediction of 
iMihos & Hernquisi (Il994ah with our improved numerical 
models: namely, that tidal torques in major mergers of gas- 
rich disks channel gas into the central regions of the galaxy, 
where it forms a dense central starburst. The starburst leaves 
a central light excess with a high phase space density, mak- 
ing the remnant more compact, reconciling the physical and 
phase space densities of disks and ellipti cal galaxies, and im- 
printin g radial gradients into the remnant (iMihos & Hernquisa 
Il994bh . 

Stars in the remnant can be separated into three distinct pop- 
ulations. First, those stars which form in the progenitor disks 
before the final merger and coalescence of the two galaxies. 
The final merger scatters the orbits of these stars as they ex- 
perience violent relaxation. They dominate the light, even in 
extremely gas-rich merger remnants, outside of 0.5 - 1 kpc, 
and develop a Sersic-law n, 2.6 ± 0.7 profile owing to their 
partial violent relaxation. These stars are effectively the dis- 
sipationless component of the merger, and behave to lowest 
order much as they would entirely in the absence of gas. For 
example, the effective radius of just these stars reflects that 
of the progenitor disks (it is substantially larger than that of a 
typical elliptical of the same mass). 

Second is the starburst population, formed in the central 
gas concentration during the final stages of the merger. Since 
dissipation allows the gas to lose energy rapidly, this compo- 
nent is very compact, and dominates the light inside a small 
radius < 0.5-1 kpc, comparable to the observed scales of cen- 
tral starbursts in merging LIRGs and ULIRGs. These stars do 
not undergo violent relaxation - because the dynamical times 
in the center of the galaxy are short, the final coalescence has 
largely completed in these central regions by the time the star- 
burst nears maximum. In other words, these stars form in a 
compact, dissipational starburst in a nearly fixed background 
potential set by the dissipationless component of the merger 
We show that the size of this dissipational component is set 
primarily by the radius at which it becomes self-gravitating. 
The gas is then generally stable against further collapse (even 
with cooling) and rapidly forms stars. Since the star formation 
time scale and subsequent collapse are both regulated by the 
dynamical time (and the system is self-gravitating), the cen- 
tral component systematically processes its gas as it contracts. 
It's worth noting that merger-induced starbursts may not be 
the only source of dissipation (for example, stellar mass loss 
may replenish the gas supply and lead to new dissipational 
bursts, see e.g. Ciotti & Ostriker 2007), but for our purposes, 
all dissipational star formation will appear similar when ob- 
served and have the same effects (we are essentially measur- 
ing the integrated amount of dissipation). 

Third, some gas survives the merger This is primarily ma- 
terial which has been moved to large radii temporarily, ei- 
ther blown out by a combination of supernova and AGN feed- 
back or thrown out in tidal tails. This gas will then slowly 
settle back in, against the background of a largely relaxed 
remnant. Since there are no strong tidal torques remaining, 
the material quickly settles into a small, rotationally sup- 
ported disk, and typically forms embedded kinematic compo- 
nen ts (embedded disks, kinernatically decoupled co res, etc.; 
see iHernquist & BarnesI [19911 iHoffman et al.ll2007l) . While 
potentially important for the kinematics of the remnant (ro- 
tation and isophotal shapes, in particular), this component 
rarely contributes substantially to the surface density profile 
of most objects, and therefore is not of immediate interest. 
Furthermore, this component has little effect on the key pre- 



dictions here, regarding e.g. the degree of central dissipa- 
tion and central mass concentration. Some caution should 
be taken, though, since in optical bands a younger embed- 
ded disk can appear much more prominent (for example, up 
to ~ 20% of the 5-band light, despite being only ~ 2% of the 
stellar mass). 

In principle, observed merger remnant light profiles can be 
decomposed into a central concentration tracing the dissipa- 
tional starburst and an outer Sersic profile reflecting the dis- 
sipationless, violently relaxed stars. However, we show that 
any attempt to infer the extra light content of a galaxy requires 
care. Fitting the light profile of a typical, albeit very gas-rich 
merger remnant to a pure Sersic law or a cored Sersic law, for 
example, and comparing the central regions with the extrapo- 
lation of the outer profile yields physically meaningless values 
of both the outer Sersic index and the extra light fraction (in 
the cored Sersic case, the fit actually is biased towards infer- 
ring that these gas-rich merger remnants have missing light in 
their centers, i.e. some kind of central mass deficit). The val- 
ues do not at all reflect the physical values: i.e. the amount 
of mass involved in the dissipational starburst or the profile of 
the violently relaxed (non-starburst) stars. That is not to say, 
however, that these are uniformly poor fits to the profiles - if 
one wishes to use such profiles to recover e.g. the effective ra- 
dius or total light content of the galaxy, the bias is not severe 
(but see Boylan-Kolchin et al. 2005). However, it highlights 
the fact that any parameterized fit will have degeneracies and 
can recover systematically different values even for a nomi- 
nally similar parameter (such as the outer Sersic index). 

By using our simulations as a testing ground, we can cal- 
ibrate such a fit in order to design a parametric decomposi- 
tion of the profile that recovers the physical values of interest. 
We find that a simple means of doing so is to fit the surface 
brightness profiles to the sum of an inner exponential (Sersic 
n = 1) and outer Sersic (free n) profile. The choice of n = 1 
for the inner regions is not meant to say that this reflects the 
true shape of the central extra light extrapolated inwards to 
very small radii, which can b e complex (with e .g. stellar clus- 
ters and features at ^ 50pc iLauer et al.|[2005h . In fact, this 
fitting procedure should not be used or extrapolated to within 
~ 30-50pc, which neither our simulations nor the observa- 
tions to which we compare typically resolve (indeed, in ex- 
periments with some well-resolved profiles, we find that in- 
cluding these small radii when fitting can lead to misleading 
results and a much higher rate of catastrophic failures of the fit 
owing to the presence of extremely small-scale features that 
are unimportant for the overall profile). The reason for the 
choice of n = 1 for the inner component is that it, on average, 
yields the correct decomposition in both total mass and radius 
between the starburst and non-starburst stars, and because it 
minimizes the degeneracy with the outer Sersic index (fitting 
an inner n = 4, for example, would introduce a large degener- 
acy with an outer profile that had a similar n ~ 4 profile). 

We apply this decomposition to our simulations and to a 
large sample of gas-rich merger remnants observed by RJ04, 
ranging from ULIRGs to shell ellipticals. We also directly 
fit each of the observed profiles to a suite of simulations - 
i.e. determine the simulation mass profile which most closely 
matches that observed. We find that, in all cases, we have 
simulations which provide good matches to the observed sys- 
tems, to better than the typical point-to-point variance inher- 
ent in the simulation surface brightness profiles (A/x < 0.1). 
We also find that the physical starburst components in these 
best-fitting simulations are closely related to those that we fit 
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directly to the observed profiles, lending further support to our 
procedure for decomposing the profiles. 

Given our analysis, we can then study how the properties 
of the two components scale. For the outer profiles, we find 
that their Sersic indices are remarkably constant as a function 
of stellar mass or any other properties. Indeed, there appears 
to be a typical Sersic index 2.5 - 3 for gas-rich merger rem- 
nants, with a scatter of ^ 0.7 about this median. We empha- 
size that this outer Sersic index is only meaningful in the sense 
of reflecting those violently relaxed stellar populations recov- 
ered in our two component decomposition. Fitting the entire 
profile to a single Sersic index can yield a very different result, 
and can introduce systematic trends (if e.g. the typical extra 
light fraction or size changes with mass). Given our attempt to 
carefully separate these components, then, this should not be 
surprising: the dissipationless component is simply acting un- 
der the influence of gravity, a nd is therefore comple tely self- 
similar across scales (see also lNaab & TruiiUoll2006h . Never- 
theless, we make a prediction that gas-rich merger remnants 
should have this approximately constant Sersic index distribu- 
tion. Subsequent gas-poor remergers may have different Ser- 
sic indices, as might pseudobulges or other low-mass bulges 
(and together these may drive a systematic dependence of Ser- 
sic index on mass, owing to the cosmological dependence of 
formation history on mass); but this particular class of ellipti- 
cals should have roughly a fixed Sersic index distribution. 

The properties of the extra light component also scale in 
a regular fashion. For a given initial simulation gas fraction, 
the extra light content of the remnant systematically decreases 
with mass. This is expected, since the star formation effi- 
ciency is higher in larger mass systems, so they consume more 
of whatever gas they have at earlier stages (e.g. on first pas- 
sage), before the final merger. What the extra light mass re- 
flects is the gas content available and channeled to the center 
during the final coalescence of the two nuclei. Comparing 
the estimated extra light masses from our two-component de- 
composition to this physical value, we find a good correlation. 
There is substantial scatter - for a given object, the inferred 
extra light fraction based on fitting the profile can be mislead- 
ing (in the sense of not reflecting the true starburst mass frac- 
tion) by a factor ^ a couple, but, on average, the appropriate 
value is recovered. 

In a few observed systems, detailed stellar population stud- 
ies have enabled estimates of how much mass was recently 
formed in a merger-induced starburst (as opposed to more ex- 
tended star formation prior to the merger). For these cases, 
then, we can directly compare our estimated extra light com- 
ponent masses to the young stellar populations, and find that 
indeed they trace one another to within a factor 2 - 3, com- 
parable to the expected scatter in both estimators. Again, this 
suggests that we can in fact infer physical decompositions 
from observed merger remnant profiles. 

There are a number of interesting applications of this de- 
composition, beyond noting the importance and physical na- 
ture of extr a light in recent gas-rich merger remnants. In 
iHopkins eT al. (2008b) we apply this analysis to large samples 
of old, relaxed elliptical galaxies with central cusps, and con- 
sider the role of the extra light in shaping the global kinematic 
properties of the galaxies. But we show here that one can 
in fact use this light component to infer something physical 
about the formation of the galaxy. Furthermore, we demon- 
strate in detail how the different components relevant to the 
surface brightness profiles of gas-rich merger remnants form 
and evolve. We also note that, with this physical decomposi- 



tion in place, we often infer extra light fractions/masses and 
radii which are much larger than previous estimates (i.e. that 
extra light is typically ^ 1 -5% of the galaxy light and be- 
comes important only within ~ Q.QSRe). This is because the 
extra light component, as it becomes larger, typically blends 
in more smoothly with the outer profile, and does not neces- 
sarily appear as a sharp departure from the outer light pro- 
file. We see this in both simulations and observations - but 
in almost every such case there is additional evidence that the 
transition to extra light is real, including changes in the stellar 
populations, ellipticity, and boxy or diskiness of the remnants 
at these radii. 

We find that gas-rich merger remnants do have excess light 
properties similar to those in gas-rich merger simulations. 
Given a careful two-component decomposition of the surface 
brightness profile (as opposed to fitting the entire profile to a 
single Sersic index, for example), we identify a statistically 
significant extra light component in every gas rich merger 
remnant observed (and simulated), with mass fractions span- 
ning a wide range ^ 3-30%. There are good simulation ana- 
logues to each observed merger remnant, and they similarly 
contain this range of masses involved in their final merger- 
induced starbursts. These extra light masses correspond to 
reasonable, expected initial gas fractions for the merging disks 
(although there is no one-to-one correspondence ). They are 
also co mparable to the estimated fractions in Hern quist et al.l 
(1 19931) . namely the dissipational mass fraction needed to ex- 
plain the discrepancies between the maximum phase space 
densities and surface brightness of ellipticals and their pro- 
genitor spiral galaxies. We therefore confirm in observations 
of gas-rich merger remnants the long-standing theoretical pre- 
diction that sufficient dissipation is required to explain this 
discrepancy in the central profiles of ellipticals. 

We have studied these properties and identified robust 
trends across a large suite of simulations, in which we vary 
e.g. the galaxy masses, initial gas fractions, concentrations, 
halo masses, presence or absence of bulges, presence or ab- 
sence of black holes, feedback parameters from supernovae 
and stellar winds, orbital parameters and disk inclinations, 
and mass ratios of the merging galaxies. This range of param- 
eters allows us to identify the most important physics. Most 
of these choices, for example, affect the surface brightness 
profile, extra light mass and radius of the extra light, concen- 
tration and effective radius of the remnant, and even its ellip- 
ticity and isophotal shape only indirectly. Ultimately, what 
determines the structure of the remnant (insofar as the proper- 
ties we have considered) is, to first order, how much mass 
is in the dissipationless (violently relaxed) component ver- 
sus the dissipational/starburst component at the time of the 
final coalescence of the merging galaxies. Therefore, varying 
e.g. the orbital parameters or initial galaxy structure can alter 
the remnant substantially, but predominantly only insofar as it 
changes the amount of gas which will be available at the time 
of the final coalescence of the galaxy nuclei (i.e. how much 
mass ends up in the starburst component, as opposed to being 
violently relaxed in the merger). 

In the simplest possible scenario outlined above, the cen- 
tral light and dissipationless component arise simultaneously 
in a single gas-rich merger. In practice, the situation need not 
be so simple. It is well-established that spheroids undergo- 
ing subsequent mergers will conserve rank order in particle 
(stellar) binding energy and (correspondingly) mass profile 
shape (Barnes 1992). As a consequence, we expect that re- 
mergers or "dry" mergers will conserve the central light ex- 
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cesses (starburst component) originally established in a gas- 
rich merger, even if the entire galaxy profile expands by a 
typical factor ~ 2 as it doubles in mass. In 'H opkins et alj 
( l2008ab . we explicitly confirm in numerical simulations that 
the distinction between the remnant excess light and dissipa- 
tionless outer profile is conserved in successive re-mergers, 
and that our methodology continues to reliably separate the 
original dissipational (compact, merger-driven starburst) and 
dissipationless (disk) components. This applies just as well 
to bulges in early-type galaxies: regardless of their formation 
mechanism (whether formed "in situ" by disk instabilities or 
minor mergers, or formed by earlier major mergers with a re- 
accreted disk), they will fundamentally be composed of a dis- 
sipationless (scattered) component and a dissipational (star- 
burst driven by energy and angular momentum loss in gas) 
component, and these components will be separately con- 
served in subsequent mergers. For any merger remnant or el- 
liptical galaxy, therefore, the observed central or excess light 
component should most generally be thought of as the sum 
of all compact starburst components formed in the history of 
the galaxy, and the outer light likewise as the sum of dissi- 
pationlessly scattered disk stars. In other words, if a spiral is 
transformed to a progressively earlier type (and, eventually, a 
true elliptical) by, say, a rapid series of many 5 : 1 mergers, 
each of which disrupts some disk stars (scattering them into a 
spheroid) and causes some gas to lose angular momentum and 
produce a small nuclear starburst (building up the central light 
component), then the final inner and outer light components 
will reflect the sum of the dissipational and dissipationless 
events. 

Our physical interpretation of the extra light, and our con- 
clusions and comparisons in this work, are therefore not 
changed by the merger histories of observed galaxies, nor by 
e.g. whether or not the progenitor disks have bulges - it is sim- 
ply possible that some of the central mass concentration was 
built up in multiple previous events along with the most recent 
merger (but the total gas/dissipational content involved is con- 
served). Of course, in modeling those progenitor bulges, this 
raises the important question of how their structure should be 
initialized: some choice must be adopted (just as if one were 
to construct an "initial" elliptical to model a "dry" merger) 
for the initial mass fraction in a very compact dissipational 
component (which will, then, remain compact or "extra" light 
in the final merger remnant) versus the more extended dissi- 
pationless component (which will become part of the dissi- 
pationless component of the final remnant). Fortunately, to 
the extent that (at least "classical") bulges obey the same pa- 
rameter correlations as ellipticals, this suggests that we can 
estimate such properties and what they imply for the original 
progenitors and their gas content by restricting our study to 
bulge-dominated systems or pure ellipticals. 

We note that there is considerable room for progress in 
modeling the extra light component itself - including its kine- 
matics, the shape of the surface brightness profile at small 
radii (<C lOOpc), and the structure of galactic nuclei near a 
central black hole. Observations are rapidly making progress 
in this area, and revealing new insights into the formation his- 
tories of elliptical galaxies. Unfortunately, modeling these 
radii in a meaningful sense will require fundamental improve- 
ments in numerical simulations. At present, our spatial reso- 
lution reaches ~ 30 pc. In principle, it would not be hard to 
improve this to < lOpc. However, at these scales, we ap- 
proach the sizes of structure in the ISM - i.e. individual giant 
molecular clouds, star forming regions, and massive star clus- 



ters or galactic stellar nuclei. Without simulations which can 
self-consistently form these structures (i.e. include the multi- 
ple gas phases of the ISM and their exchange), as well as re- 
solve e.g. individual supernova blastwaves and remnants, im- 
proved spatial resolution has no physical meaning. It is there- 
fore an important and ambitious goal that the next generation 
of studies of galactic nuclei move beyond the sub-resolution 
prescriptions necessary when modeling large scales and at- 
tempt to include star formation, supernova feedback, and re- 
alistic, resolved ISM structure in simulations of galaxy merg- 
ers. 

In eariier work (e.g. iHopkins et all l2006alfbL l2008clldh we 
developed a model linking starbursts, quasar activity, the 
growth of supermassive black holes, and the origin of el- 
lipticals through evolutionary phases of the same events, 
driven by mergers of gas-rich galaxies. There is much 
observational support for the various links in this chain. 
ULIRGs are invariably a ssociated with gas-rich mergers (e.g. 
ISanders & MirabellI996l) and have bolo metric pro perties sim- 
ilar to bright quasars (e.g. lSanders et alJ ll988a b"). suggesting 
that ULIRGs evolve into q uasars. By the Soltan argument 
(ISoltanlll982t iHopkins et a l. 2007d), the bulk of the cosmic 
mass density in supermassive black holes was accumulated 
during periods of bright quasar activity. Observed correlations 
between supermassive black holes and properties of their host 
elhpticals (Magorrian et al. 1998; Ferrarese & Merritt 200(1 
iGebhardlTet al. 2000) demonstrate that they formed together, 
not independently. To the extent that gas-rich mergers were 
responsible for growing most of the mass in supermassive 
black holes, these correlations strongly endorse the view that 
ellipticals formed originally in gas-rich mergers of spirals. 
(With the possibility that they could have been modified sub- 
sequently in gas-free mergers with other ellipticals that did not 
trigger quasar activity or lead to black hole growth, and that 
modest black hole growth can occur with non-merger induced 
fueling (Hopkins & Hernquist 2006).) 

In addition to these observational lines of evidence, there 
are simple physical arguments that support these connections. 
The bulk of the stellar mass in ellipticals was l ikely assem- 
bled d issipationlessly through violent relaxation (iLynden-Belll 
Il967h . but the growth of supermassive bla ck hol es proba- 
bly involved accretion of gas (lLvnden-Beliill969h . As we 
argue m e.g. I Hopkins et all (l2008dh . the condition that su- 
permassive black holes and ellipticals originate together re- 
quires their assembly in gas-rich mergers, which contain a 
galaxy's worth supply of both gas and stars. Our analysis 
here lends additional credence to this hypothesis: the outer 
components of the light profiles of ellipticals were indeed put 
into place through self-similar gravitational physics, but the 
compact, inner starburst populations resulted from gas dissi- 
pation. The same gas dissipation that yielded this inner com- 
ponent also provided the material to grow supermassive black 
holes in ellipticals in a self-regulated manner, accounting for 
the similarities between the fundamental plane of ellipticals 
(iDressler et al.ll987HDjorg ovski & Davis"198'A and the black 
hole fundamental plane (Hopkins et al. 2007a b). Observa- 
tionally, then, we see that this blending of dissipationless stel- 
lar dynamics and gas dissipation is reflected not only in corre- 
lations between supermassive black holes and their hosts, but 
in the detailed structure of elliptical galaxies as well. 
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APPENDIX 
FITS TO THE SAMPLE OF RJ04 

In Figures I2TII38] we explicitly show the results of our fitting and simulation comparison with each of the merger remnants in 
the sample of RJ04. Table[T]summarizes the results, including the estimated integrated properties from simulations corresponding 
to each observed system. 
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Fig. 21. — A.S Figure|6] Figures l2U38l show the results for all objects in the RJ04 sample. Objects are sorted from most to least luminous in A'-band. Note 
that Mrkl014 (previously unpublished) and UGC8058 (Mrk231) shown here are contaminated by central AGN, giving rise to the large discrepancies and poor 
fits seen. 




Fig. 22. — Figure [2n continued. 
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Fig. 23. — Figure [2T] continued. 




Fig. 24. — Figure [2n continued. 
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Fig. 25. — Figure llT] continued. 
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Fig. 26. — Figure [2n continued. Unrelaxed features (including remaining spiral and tidal structure) in UGC9829 prevent a reliable fit. 
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Fig. 27. — Figure lzT] continued. 




Fig. 28. — Figure [2n continued. 



Extra Central Light in Merger Remnants 3 1 




Fig. 29. — Figure llT] continued. 
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Fig. 3 1 . — Figure [2T] continued. 




Fig. 32. — Figure llTI continued. 
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Fig. 33. — Figure llT] continued. 




Fig. 34. — Figure llTI continued. 
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Fig. 35. — Figure [2T] continued. NGC6052 appears to be still unrelaxed at ~a few kpc. 




Fig. 36. — Figure [211 continued. There is a thin, bar-like feature at ~ 1 — 2kpc in NGC4004, not reproduced in our simulations at the observed snapshots 
(although similar features do generically occur in earlier merger stages). 
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Fig. 37. — Figure llT] continued. NGC3310 shows a face-on ring and arm structure at ~ 1 kpc, pushing the simulated extra light component to larger radii. 




Fig. 38.— 



Figure lzTI continued. 
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TABLE 1 
Fits to RJ04 Merger Remnants 



Name 


R.A. (J2000) 


Decl. (J2000) 


Mk 


n, (fit)' 


Hj (sim)^ 


n, (RJ04)3 


fe (fit/ 


fsb (sim)' 


/gas (sim)* 


A/i (fit) 


A/i (sim; 


Mrkl014^ 


01 59 50 


00 23 41 


-28.16 


3.9 


4.7 ±3.2 


10.0 


0.67 


0.24±0.15 


0.40 ±0.21 


0.06 


0.24 


UGC8058'' 


12 56 14 


56 52 25 


-27.55 


4.6 


1.6 ±0.9 


10.0 


0.11 


0.37 ±0.14 


0.49 ±0.25 


0.14 


0.42 


Arpl56 


10 42 38 


77 29 41 


-25.81 


3.1 


3.0 ±1.0 


9.64 


0.12 


0.17 ±0.06 


0.24 ±0.19 


0.11 


0.06 


AM06 12-373 


06 13 47 


-37 40 37 


-25.65 


1.6 


2.5 ±0.9 


3.44 


0.16 


0.08 ± 0.04 


0.05 ± 0.05 


0.07 


0.06 


AM2246-490 


22 49 39 


-48 50 58 


-25.52 


2.8 


3.1 ±2.0 


10.0 


0.32 


0.19 ±0.07 


0.49 ±0.16 


0.12 


0.10 


NGC6598 


18 08 56 


69 04 04 


-25.51 


3.4 


3.2± 1.1 


3.65 


0.05 


0.10 ±0.06 


0.05 ± 0.06 


0.05 


0.04 


UGC5101 


09 35 5 1 


61 21 11 


-25.50 


4.6 


2.4 ±0.7 


10.0 


0.32 


0.45 ±0.14 


0.54 ±0.19 


0.04 


0.04 


NGC828 


02 10 09 


39 11 25 


-25.36 


3.6 


3.2 ±0.7 


2.99 


0.19 


0.24 ±0.16 


0.40 ± 0.23 


0.07 


0.08 


NGC2418 


07 36 37 


17 53 02 


-25.31 


2.1 


2.6 ±0.8 


2.85 


0.13 


0. 1 1 ± 0. 1 1 


0.16 ±0.23 


0.05 


0.09 


Arpl87 


05 04 53 


-10 14 51 


-25.25 


3.9 


2.5± 1.1 


4.10 


0.06 


0.12 ±0.08 


0.14 ±0.21 


0.10 


0.08 


UGC10607 


16 55 09 


26 39 46 


-25.20 


1.7 


2.9± 1.2 


2.94 


0.33 


0.30 ±0.17 


0.48 ± 0.20 


0.10 


0.09 


NGC5018 


13 13 00 


-19 31 05 


-25.15 


3.1 


3.7 ±2.0 


4.39 


0.06 


0.10 ±0.06 


0.34 ±0.20 


0.07 


0.07 


NGC3921 


11 51 06 


55 04 43 


-25.13 


2.5 


2.6 ±0.4 


5.03 


0.26 


0.27 ±0.13 


0.50±0.15 


0.09 


0.06 


AM03 18-230 


03 20 40 


-22 55 53 


-25.09 


2.3 


4.6 ±2.0 


5.01 


0.28 


0.25 ±0.11 


0.49 ±0.14 


0.13 


0.06 


AM2055-425 


20 58 26 


-42 39 00 


-25.08 


0.8 


3.3 ±3.5 


2.34 


0.34 


0.25 ±0.13 


0.42 ± 0.20 


0.16 


0.1 1 


NGC7585 


23 18 01 


-04 39 01 


-24.98 


2.4 


3.3±0.8 


3.53 


0.11 


0.10 ±0.06 


0.18 ±0.20 


0.05 


0.06 


UGC9829^ 


15 23 01 


-01 20 50 


-24.96 


1.4 


3.2±1.1 


2.34 


0.02 


0.09 ± 0.08 


0.13 ±0.20 


0.32 


0.22 


AM1419-263 


14 22 06 


-26 51 27 


-24.94 


3.2 


3.7 ±3.5 


4.12 


0.06 


0.10 ±0.03 


0.12±0.18 


0.05 


0.05 


AM1255-430 


12 58 08 


-43 19 47 


-24.93 


0.9 


4.7 ±2.2 


1.87 


0.26 


0.46 ±0.11 


0.40 ± 0.24 


0.09 


0.07 


IC5298 


23 16 00 


25 33 24 


-24.92 


1.6 


2.6 ±3.5 


3.93 


0.34 


0.30±0.12 


0.49 ±0.18 


0.14 


0.14 


NGC7252 


22 20 44 


A Ai\ A -\ 

-24 40 41 


-24.84 


1.3 


3.0 ±0.5 


3.32 


0.32 


0.23 ±0.14 


0.49 ±0.21 


0.06 


0.09 


UGC10675 


17 03 15 


31 27 29 


-24.80 


2.2 


4.6 ±2.8 


10.0 


0.64 


0.51 ±0.08 


0.53±0.18 


0.06 


0.10 


NGC1614 


04 33 59 


-08 34 44 


-24.74 


2.3 


3.1 ±0.9 


10.0 


0.35 


0.34 ±0.11 


0.50 ±0.12 


0.14 


0.14 


NGC3256 


10 27 51 


-43 54 14 


-24.72 


1.8 


2.9 ±0.8 


2.17 


0.09 


0.35 ±0.12 


0.49 ±0.24 


0.11 


0.15 


UGC4635 


08 51 54 


40 50 09 


-24.71 


2.9 


2.7 ±0.5 


5.07 


0.08 


0.11 ±0.05 


0.49 ±0.21 


0.22 


0.13 


AM2038-382 


20 41 13 


-38 11 36 


-24.70 


2.5 


2.4 ± 1.1 


10.0 


0.50 


0.43 ± 0.08 


0.51 ±0.15 


0.08 


0.07 


AM 1300-233 


13 02 52 


-23 55 18 


-24.65 


1.9 


3.4 ±2.0 


4.99 


0.10 


0.10 ±0.08 


0.12±0.18 


0.19 


0.17 


NGC455 


01 15 57 


05 10 43 


-24.64 


2.0 


2.7 ± 1.1 


6.21 


0.20 


0.19 ±0.04 


0.34 ± 0.20 


0.06 


0.06 


NGC34 


00 1 1 06 


-12 06 26 


-24.61 


3.5 


2.4 ± 1.4 


10.0 


0.39 


0.51 ±0.09 


0.51 ±0.17 


0.12 


0.06 


UGC2238^ 


02 46 17 


13 05 44 


-24.58 


1.1 


4.7 ± 2.0 


1.46 


0.16 


0.36 ±0.11 


0.47 ± 0.24 


0.13 


0.20 


UGC11905 


22 05 54 


20 38 22 


-24.51 


1.2 


2.6±3.5 


5.24 


0.42 


0.30±0.10 


0.51 ±0.15 


0.06 


0.09 


Arpl93 


13 20 35 


34 08 22 


-24.40 


1.2 


3.2±3.3 


2.74 


0.33 


0.30 ±0.11 


0.49 ±0.20 


0.09 


0.15 


NGC7727 


23 39 53 


-12 17 35 


-24.23 


2.6 


3.0±0.6 


3.41 


0.13 


0.23 ±0.14 


0.48 ±0.23 


0.08 


0.08 


NGC2623 


08 38 24 


25 45 17 


-24.22 


4.9 


3.3 ± 3.4 


10.0 


0.17 


0.31 ±0.08 


0.51 ±0.09 


0.18 


0.13 


UGC6 


00 03 09 


21 57 37 


-24.01 


2.1 


2.4 ±0.7 


10.0 


0.63 


0.45 ± 0. 1 


0.51 ±0.15 


0.04 


0.11 


JNLrC / i JJ 


11 4y 4o 


—34 52 35 


—23.95 


4. 1 


3.7 ± 1.9 


10.0 


0.08 


n 1 /I J_ n HA 
U. 14 zh U.Uo 


O.Zi ± O.ZO 


0. 14 


0. 10 


JNuC27o2 


Uy 14 05 


4U Uo 4y 


—23.83 


1.7 


2.6±2.1 


6.68 


0.23 


u.iy ± U.U7 


U.37 ± u.iy 


0.13 


0.12 




f\l f\A. 

xjl DD UO 


JD 4z 13 


—Id. Io 


9 1 


Z.O ± l.U 


dMo 


U.U4 


U.Uo iL U.U4 


U.UD iL U.U/ 


u. tu 


u.uy 




11 14 41 






1^ 
U.O 


2 4± 1 4 


1 R7 


L/.Dl 




n ^fi 4- n 1 R 




U. ID 




CV\ Ofi A'\ 

\JJ WD 


—9^ AO '^0 




1 Q 




A OR 




n 1-2 4-0 14 


n AQ 4- n 99 


no 


U.UD 




1 1 XJ> JO 


SO "^0 

Oj J\j JyJ 


— 9'^ 70 


A 
o.f 




4 47 


OA 


n 1 a 4- fi OR 


n "^d. 4- n 90 




12 


NGC2655 


08 55 37 


78 13 23 


-23.70 


2.4 


3.2 ±0.9 


3.01 


0.05 


0.41 ±0.13 


0.50 ±0.17 


0.04 


0.09 


NGC60528 


16 05 12 


20 32 32 


-23.55 


0.9 


1.7 ±2.9 


1.13 


0.04 


0.10±0.15 


0.13 ±0.22 


0.11 


0.21 


NGC2914 


09 34 02 


10 06 31 


-23.51 


0.6 


5.1 ±3.4 


6.97 


0.50 


0.25 ±0.09 


0.49 ±0.13 


0.11 


0.14 


NGC4194 


12 14 09 


54 31 36 


-23.21 


1.6 


2.4 ±0.5 


4.59 


0.53 


0.47 ±0.10 


0.51 ±0.15 


0.14 


0.16 


NGC4441 


12 27 20 


64 48 06 


-22.98 


2.5 


2.6 ±2.6 


6.83 


0.14 


0.12±0.06 


0.36±0.19 


0.05 


0.06 


NGC40048 


1158 05 


27 52 44 


-22.89 


1.5 


2.6 ±0.8 


1.53 


0.42 


0.23 ±0.11 


0.32 ±0.24 


0.25 


0.19 


NGC2744 


09 04 38 


18 27 37 


-22.83 


1.7 


3.2 ±0.7 


2.35 


0.05 


0.19±0.16 


0.12 ±0.22 


0.06 


0.08 


AM1158-333 


12 01 20 


-33 52 36 


-22.61 


3.0 


3.1 ±1.0 


3.75 


0.32 


0.31 ±0.09 


0.50 ±0.13 


0.09 


0.08 


NGC33108 


10 38 45 


53 30 05 


-22.07 


1.6 


2.4±1.0 


2.41 


0.08 


0.35 ±0.12 


0.49 ±0.14 


0.20 


0.23 


Arp230 


00 46 24 


-13 26 32 


-21.75 


1.2 


2.7 ±0.6 


1.56 


0.06 


0.43 ±0.15 


0.50±0.17 


0.17 


0.13 


AM0956-282 


09 58 46 


-28 37 19 


-20.50 


2.2 


2.6 ±0.8 


2.45 


0.21 


0.23 ±0.06 


0.47 ±0.14 


0.03 


0.06 



Outer Sersic index its of the two-component best-fit profile. Range of outer Sersic indices fit in the same manner to the best-fit simulations, at / 1 - 3 Gyr after the merger when 
the system has relaxed.^ Sersic index fit to the entire profile (i.e. not decomposed into an outer Sersic and inner extra light component) in RJ04. Note the authors impose a maximum 
n,i = 10.0.'* Fraction of light in the inner or "extra light" component of the fits.^ Fraction of light from stars produced in the central, merger-induced starburst in the best-fit simulations 
(± the approximate interquartile range allowed). Note that this may be biased to high values for actively star-forming (especially LIRG and ULIRG) systems.^ Initial gas fraction in the 
best-fitting simulations - this is less robust (note the large interquartile range), representing a rough gas fraction of the systems ~a few Gyr before the final merger, if they evolved in 
isolation. AGN contamination affects the central regions (i.e. estimated extra light) and prevents a good simulation fit (we only model the stellar profile here). ^ Unrelaxed or prominent 
disk/bar features make the comparison with these objects uncertain. 



